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ABSTRACT 


Naval  power  distribution  has  principally  used  an  AC  network  to  supply  loads. 
With  the  advent  of  new  power  electronic  devices,  the  focus  has  shifted  to  employing  a 
DC  distribution  system  that  eliminates  large  transformers  and  mechanical  switching 
devices  and  enhances  the  survivability  of  the  platform.  The  Power  Electronic  Building 
Block  (PEBB)  Network  Simulation  Testbed  currently  under  construction  at  the  Naval 
Postgraduate  School  is  a  study  into  the  feasibility  of  such  DC  systems. 

The  objective  of  this  thesis  was  to  perform  theoretical  and  simulation-based 
analysis  to  establish  quantitative  criteria  for  PEBB  Testbed  stability.  These  criteria  were 
then  used  to  develop  a  set  of  hardware  studies  to  investigate  the  interaction  of 
components  within  the  PEBB  Testbed.  Finally,  the  hardware  studies  were  utilized  to 
verify  PEBB  Testbed  performance. 

Principal  conclusions  of  this  research  included  that  the  PEBB  Testbed 
demonstrated  stability  under  all  simulated  loading  conditions.  Follow-on  testing  of  the 
PEBB  Testbed  confirmed  that  the  simulations  correlated  well  with  hardware 
implementation.  In  addition,  the  hardware  validation  studies  revealed  that  switching 
harmonics  had  a  considerable  effect  on  the  system  output. 
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I.  INTRODUCTION 

Naval  power  distribution  has  principally  used  AC  networks  to  supply  loads.  With 
the  advent  of  new  power  electronic  devices,  the  focus  has  shifted  to  employing  a  DC 
distribution  system  that  eliminates  large  transformers  and  mechanical  switching  devices 
and  enhances  the  survivability  of  the  platform.  The  Power  Electronic  Building  Block 
(PEBB)  Network  Simulation  Testbed  currently  under  construction  at  the  Naval 
Postgraduate  School  will  be  used  to  study  the  feasibility  of  such  DC  systems. 

The  proposed  architecture  for  this  shipboard  power  distribution  scheme  is  shown 
in  Figure  (1-1).  In  this  distribution  network  the  main  feeders  are  DC.  The  ship  is 


Figure  1-1,  Integrated  Power  System  (from  [Ref.  1]) 
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divided  into  zones  with  each  zone  containing  common  energy  conversion  devices  fed 
from  the  DC  busses.  The  DC  power  is  distributed  via  port  and  starboard  busses  from  the 
source(s)  into  the  separate  zones.  Each  zone  contains  a  number  of  Ship  Service 
Converter  (SSCM)  and  Inverter  Modules  (SSIM).  The  SSCM  is  used  in  each  zone  to 
step-down  the  distribution  bus  voltage  to  a  regulated  level  for  use  in  the  zone.  In  this 
way  the  SSCM  inserts  intelligence  into  the  system  by  acting  to  buffer,  preregulate  and 
fault  protect  each  zone.  Electric  loads  within  the  zone  are  fed  by  either  SSCMs  or  SSIMs 
depending  on  the  load’s  requirement  for  DC  or  AC  power. 

The  focus  of  this  thesis  is  on  the  design  of  SSCMs  and  their  subsequent 
integration  with  SSIMs  in  a  DC  distribution  network.  Coincident  to  this  research  is  the 
development  of  system  stability  criteria.  As  a  result,  theoretical  and  simulation-based 
analysis  will  be  performed  to  establish  quantitative  criteria  for  system  stability.  These 
criteria  will  then  be  used  to  develop  a  set  of  hardware  studies  to  investigate  the 
interaction  of  components  within  the  PEBB  Network  Simulation  Testbed.  Finally,  the 
hardware  studies  will  be  conducted  to  verify  system  performance. 

The  documentation  detailing  this  research  is  organized  into  six  chapters.  In 
Chapter  II,  the  basic  design  of  the  SSCM  power  section  is  presented.  Chapter  III  deals 
with  the  development  of  SSCM  analog  and  digital  closed-loop  state-difference  control 
algorithms.  Pole  placement  and  gain  selection  for  these  controllers  are  investigated  in 
Chapter  IV.  With  the  SSCM  gains  determined,  detailed  simulations  of  integrated 
SSCM/SSIM  operations  are  documented  in  Chapter  V.  This  chapter  also  includes 
hardware  study  results  for  some  of  the  simulation  configurations  studied. 
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Circuit  operations  and  schematics  for  the  SSCM  power  section  and  controller  are  covered 
in  Chapter  VI.  The  final  chapter  contains  a  summary  of  the  research  work,  notable 
conclusions  and  recommendations  for  future  work. 
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II.  POWER  SECTION  DESIGN 


A.  TOPOLOGY 


This  project  focused  on  the  design,  implementation  and  testing  of  high-bandwidth 
DC-to-DC  buck  converters.  Figure  (2-1)  illustrates  the  basic  circuit  topology. 


acts  as  a  DC  step-down  transformer.  The  reduction  in  voltage  is  governed  by  the  duty 
cycle  of  the  switch  according  to 

VC=DE  (2-1) 

where  D  is  the  duty  cycle.  From  Equation  (2-1),  it  is  evident  that  any  desired  output 
voltage  can  be  attained  by  controlling  the  duty  cycle. 

System  design  called  for  the  application  of  this  topology  as  part  of  the  DC  Zonal 
Electrical  Distribution  (DC  ZED)  Network  envisioned  for  the  Surface  Combatant  for  the 
twenty-first  century.  As  part  of  the  basic  research  into  DC  ZED,  a  Power  Electronic 
Building  Block  (PEBB)  Network  Simulation  Testbed  is  being  assembled  at  the  Naval 
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Postgraduate  School.  A  general  topology  of  this  network  is  pictured  in  Figure  (2-2). 
From  the  diagram,  it  can  be  seen  that  two  Source  Buck  Converters  condition  rectified 
three-phase  power  and  provide  the  network  with  regulated  DC  power.  This  power  is 
directed  to  loads  supplied  by  two  Load  Buck  Converters  and  four  Auxiliary  Resonant 
Commutated  Pole  (ARCP)  Inverters. 


3  Phase  AC 
Power  Supply 


ARCP#1 


Terminal 

Block 


Source 

Buck#! 


ARCP  #2 


AC  Circuit 


3  Phase 
Rectifier 


Source 
Buck  #2 


ARCP  #3 


ARCP  #4 


Load  Bode  HI 


Load  Buck  #2 


Load#  I 

Load  #2 

Load  #3 

1 

Load  M 

-  » 

Load  #5 

_ 1 _ 

Load  #6 

Figure  2-2,  DC  Distribution  System 

Specifications  for  the  Source  and  Load  Buck  Converters  are  given  below.  Detailed 
design  descriptions  of  the  ARCP  units  are  given  in  reference  [3]. 


B.  SPECIFICATIONS  AND  COMPONENT  SELECTION 

Minimum  design  specifications  for  the  Source  and  Load  Buck  Converters  were 
provided  by  the  Naval  Surface  Warfare  Center,  Carderock  Division/Annapolis 
Detachment.  Table  (2-1)  summarizes  the  requirements.  Using  this  information, 
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component  sizing  and  selection  was  performed.  The  final  results  of  this  analysis  are 
documented  in  Table  (2-2).  Analysis  of  each  component  selection  is  given  below. 


Parameter 

Source  Buck 

Load  Buck 

Maximum  Output  Power 

9  kW 

3  kW 

Switching  Frequency 

5-25  kHz  Hard  Switched 

5-25  kHz  Hard  Switched 

Input  (nominal) 

400  VDC  /  22.5  A 

300  Vdc/  10  A 

Output  (nominal) 

300  VDC  /  30  A 

208  V^/ 14.5  A 

Continuous  Operations 

10%  to  100%  Loading 

10%  to  100%  Loading 

Table  2-1,  Buck  Converter  Specifications 


Component/Parameter 

Source  Buck 

Load  Buck 

Switch  /  Diode 

International  Rectifier 
IRGT1090U06  600V/90A 
IGBT  Power  Modules 

International  Rectifier 
IRGT1090U06  600V/90A 
IGBT  Power  Modules 

Input  Filter  Capacitor 

Sprague  Powerlytic 

2000  pF,  450  Vdc 

Two  Sprague  Powerlytic 

230  pF,  450  VDC,  Paralleled 

Input  Filter  Inductor 

0.425  mH,  hand  wound 

0.425  mH,  hand  wound 

Output  Capacitor 

Sprague  Powerlytic 

400  pF,  450  VDC 

Sprague  Powerlytic 

400  pF,  450  VDC 

Output  Inductor 

0.760  mH,  hand  wound 

0.875  mH,  hand  wound 

1 .3  mH,  hand  wound 

Switching  Frequency 

20  kHz 

20  kHz 

Output  Voltage  Ripple 

<1  % 

<1  % 

Table  2-2,  Component  Selection  Results 


1.  Switch  and  Diode  Selection 

For  both  the  Source  and  Load  Buck  Converters,  International  Rectifier  600V/90A 
IGBT  Power  Modules  were  chosen  based  on  their  high  current  density,  rugged  design, 
and  simple  gate  drive.  In  addition,  the  INT-A-pak  “half-bridge”  packaging  lends  itself  to 
the  buck  converter  topology.  One  power  module  contains  two  IGBT  /  diode  pairs 
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allowing  co-location  of  the  switch  and  diode  components.  Appendix  (A)  contains  the 
data  sheet  for  this  component. 


2.  Output  Inductor  Sizing  and  Switching  Frequency 
a.  Source  Buck  Converter 

Inductor  sizing  was  done  to  maintain  continuous  inductor  current  over  a 
load  range  of  10%  ->•  100%.  Based  on  this  range,  the  most  limiting  condition  for 
continuous  inductor  current  occurs  at  minimum  load.  To  maintain  continuous  operations 

at  this  point,  the  inductor  must  be  wound  to  meet  or  exceed  the  critical  inductance  for  the 
circuit.  This  value  is  given  in  reference  [2]  by 

T .  R 

^crit  —  2  (2-2) 

where  T  is  the  switching  period  and  R  is  the  resistance  at  minimum  loading. 

Analysis  of  Equation  (2-2)  shows  that  two  of  the  input  parameters  are  set. 
From  Equation  (2-1),  the  duty  cycle  is  75%.  In  addition,  load  resistance  is  100  O  based 
on  a  nominal  output  voltage  of  300  V.  As  a  result,  only  the  switching  period,  or 
conversely  the  switching  frequency,  must  be  set.  Utilizing  Table  (2-1),  the  switching 
frequency  was  chosen  to  be  20  kHz.  This  selection  corresponds  to  a  switching  period  of 
50  psec  and  provides  the  following  advantages: 

•  minimization  of  required  inductor  size 

....  #  y 

•  minimization  of  the  inductor  ripple  current  given  by  A I  — - — .  DT 

L 

•  elimination  of  audible  switching  noise  of  the  buck  converter 
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•  maintains  a  5  kHz  margin  to  the  maximum  hard-switched  limit  of  the  IGBTs. 
Substituting  into  Equation  (2-2),  critical  inductance  was  found  to  be  625  jiH.  The  actual 
inductances  achieved  for  the  two  Source  Buck  Converters  were  760  pH  and  875  pH 
respectively. 

This  sizing  was  not  made  any  larger  than  necessary  because  of  the  tradeoff 
between  desired  inductance  and  the  size  required  to  achieve  that  inductance.  At  full  load, 
inductor  current  has  a  30  A  dc  offset.  This  dc  component  drives  a  single,  wound  inductor 
core  sized  to  L^,  into  saturation.  As  a  result,  each  inductor  was  made  from  several  cores 
wound  in  series.  Cores  were  wound  to  minimize  saturation  effects  while  maximizing 
inductance  per  core.  The  overall  inductance  is  then  given  by  summing  each  of  the  series 
core  inductances.  Reference  [4]  contains  core  sizing  information. 

b.  Load  Buck  Converter 

Inductor  sizing  was  done  to  maintain  continuous  inductor  current  over  a 
load  range  of  10%  -»  100%.  A  20  kHz  switching  frequency  was  chosen  for  the  reasons 
mentioned  above.  Substituting  appropriate  values  into  Equation  (2-2)  yielded  a  critical 
inductance  of  1.1  mH.  In  actuality,  a  1.3  mH  inductance  value  was  achieved  for  each 
Load  Buck  Converter.  Contrary  to  the  Source  Buck  Converters,  core  saturation  due  to  dc 
offset  was  not  a  concern  since  the  maximum  load  current  is  well  below  saturation  levels. 
As  a  result  each  inductor  was  wound  using  one  core  thus  saving  space  in  the  Load  Buck 
Converter  topology. 
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V 


3.  Output  Capacitor  Sizing 


a.  Source  Buck  Converter 

Capacitor  selection  was  based  on  maintaining  less  than  1%  voltage  ripple 
on  the  Source  Buck  Converter  output.  From  reference  [2],  percent  voltage  ripple  is 
given  by: 


AK 

K 


r2(  i-p) 

8  LC 


•100 


(2-3) 


Substituting  the  appropriate  values,  the  minimum  output  capacitance  is  found  to  be 

10  pF.  Based  on  this  small  size,  rapid  transient  response  can  be  expected.  Computer 

simulation  and  gain  selection  detailed  in  Chapter  IV,  however,  shows  that  the  transient 

voltage  response  is  unacceptably  soft  at  such  low  output  capacitance  values.  Utilizing 

these  results,  400  pF  Sprague  electrolytic  capacitors  were  selected.  This  choice  provided 

the  required  transient  response  and  built  in  a  design  margin  for  testing  above  the  10%  -> 

0 

100%  10%  load  transient  used  for  analysis  in  Chapter  IV. 


b.  Load  Buck  Converter 

Capacitor  selection  was  based  on  maintaining  less  than  1%  voltage  ripple 
on  the  Load  Buck  Converter  output.  Utilizing  Equation  (2-3),  minimum  output 
capacitance  was  determined  to  be  7.5  pF.  Using  computer  simulation  and  gain  selection 
results  detailed  in  Chapter  IV,  400  pF  Sprague  electrolytic  capacitors  were  selected  to 
provide  the  required  voltage  transient  response. 
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4. 


Input  Filter  Sizing 


a.  Source  Buck  Converter 

One  circuit  element  not  shown  in  the  basic  buck  converter  topology  of 
Figure  (2-1)  is  the  input  LC  lowpass  filter.  The  input  filter  is  documented  in  Figure  (6-1) 
of  Chapter  VI.  Initial  filter  design  focused  on  component  selection  to  provide  10  dB  of 
attenuation  at  60  Hz.  Component  sizing  and  dc  attenuation  concerns,  however,  make  this 
design  impractical.  As  a  result,  ac  ripple  rejection  analysis  shifted  to  a  combined 
filter/controller  loop  speed  method.  In  this  arrangement,  the  input  LC  filter  decouples  the 
rectifier  source  from  the  Source  Buck  Converter  input  by  at  least  10  dB  of  attenuation  for 
frequencies  360  Hz  and  above.  The  controller  loop  speed  of  the  Source  Buck  Converter 
is  then  tasked  with  providing  ac  ripple  rejection  below  360  Hz.  Chapter  IV  contains  a 
discussion  on  controller  loop  speed. 

Central  to  the  selection  of  capacitor  and  inductor  values  was  the  placement 
of  the  resonant  peak  of  the  input  LC  filter.  The  resonance  acts  to  amplify  low  frequency 
components  in  the  supply  voltage.  As  a  result,  peak  placement  must  minimize  this 
amplification  for  the  frequencies  of  concern.  Since  the  supply  voltage  is  rectified  3 -phase 
power,  these  frequencies  are  60  Hz,  120  Hz,  and  360  Hz.  Matlab  analysis  revealed  that  a 
resonant  peak  placed  at  180  Hz  provides  the  desired  response.  This  value  was  obtained 
using  a  capacitance  of  2000  pF  and  an  inductance  of  391  pH.  In  actuality,  a  1 73  Hz 
resonant  peak  was  achieved  using  a  2000  pF  capacitor  and  425  pH  inductor. 
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b.  Load  Buck  Converter 

The  input  filter  to  the  Load  Buck  Converters  was  designed  to  provide 
decoupling  from  the  output  of  the  Source  Buck  Converters.  As  a  result,  the  LC 
combination  had  to  be  chosen  to  place  the  resonant  peak  of  the  Load  Buck  Converter 
input  filter  below  the  slowest  controller  poles  of  the  Source  Buck  Converters.  Based  on 
the  pole  placement/gain  selections  detailed  in  Chapter  IV,  the  resonant  peak  of  the  input 
LC  filter  was  set  at  360  Hz  using  a  capacitance  of 460  pF  and  an  inductance  of  425  pH. 

C.  SUMMARY 

With  the  basic  power  sections  of  the  Source  and  Load  Buck  Converters 
established,  the  design  effort  focused  on  the  development  and  implementation  of  an 
appropriate  control  scheme.  As  will  be  shown  in  Chapter  III,  open-loop  operations  of  a 
buck  converter  are  extremely  underdamped  under  transient  conditions.  Therefore, 
satisfactory  system  transient  response  must  be  obtained  by  “closing  the  loop”.  Once  the 
required  control  algorithm  has  been  established,  gain  selections  and  simulations  can  be 
performed  to  produce  stability  criteria  for  hardware  testing. 
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III.  CONTROL  ALGORITHM  DEVELOPMENT 


A.  OPEN-LOOP  ABCD  STATE  SPACE  MODEL 

The  first  step  in  controller  development  required  the  determination  of  the  open- 
loop  state  space  equations.  In  state  space  analysis,  energy  storage  devices  dictate  the 
state  variables.  Referring  to  Figure  (2-1)  in  Chapter  II,  the  obvious  state  variable  choices 
are  capacitor  voltage  and  inductor  current.  With  the  switch  shut,  the  state  equations  are: 


_  1 
dt  C 


R 


(3-1) 


dt 


=  f[*-Vc] 


(3-2) 


With  the  switch  open,  the  state  equations  become: 


dvc 

d  t 


C 


lL 


R 


dt  Ll  CJ 


(3-3) 

(3-4) 


Combining  like  equations  into  an  averaged  model  yields 


dvc 

dt 


h~ 


R 


^  =  -[dE-vc  ] 
dt  Ll  ci 


(3-5) 

(3-6) 


where  the  duty  cycle,  d,  is  given  by: 


d  = 


^  switch  on 


t  +  t  T 

*  switch  on  switch  off  switch 


(3-7) 
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with  Tswitch  “  switching  period.  The  resulting  ABCD  state  space  representation  is: 


When  operated  in  the  open-loop  configuration,  the  buck  converter  output  voltage  and 
inductor  current  responses  are  highly  underdamped  in  the  presence  of  load  transients. 


Figure  (3-1)  illustrates  the  system  response  for  a  buck  converter  operating  at  a  75%  duty 
cycle  with  E  =  400  V,^. 


B.  KASSAKIAN  ALGORITHM 

The  large  overshoots,  long  settling  times  and  high-frequency  oscillations  shown  in 
Figure  (3-1)  are  not  acceptable  in  an  interconnected  power  system.  Clearly,  a  need  has 
been  established  to  “close  the  loop”.  The  foundation  for  the  final  control  algorithm 
comes  from  reference  [5]  and  is  given  by: 

d(t)  =  - ht  i L  (f) -hvVo(t)-hn\ v0 (g)dg  (3-10) 

This  multiloop  control  law  states  that  duty  cycle  perturbations  are  a  function  of 
perturbations  in  the  average  output  voltage  and  average  inductor  current.  By  judicious 
selection  of  the  constant  feedback  gains  hf,  hv,  and  hn,  the  desired  closed-loop  response 
can  be  obtained. 

The  last  term  in  Equation  (3-10)  represents  integral  control  action.  Assuming  that 
the  buck  converter  closed-loop  system  is  stable  and  driven  only  by  constant  signals,  all 
variables  must  settle  to  constant  values  in  the  steady  state.  As  a  result,  the  integrand  in 
Equation  (3-10)  must  settle  to  zero  to  prevent  the  integral  from  contributing  a  time- 
varying  term  to  the  equation.  -  Therefore,  so  long  as  the  feedback  gains  are  chosen  to 
stabilize  the  system,  zero  steady-state  error  is  assured  in  the  output  voltage. 

With  Equation  (3-10)  serving  as  a  model  for  the  feedback  path,  the  closed-loop 
state  equations  were  determined.  This  step  was  done  by  performing  small-signal  analysis 
of  Equation  (3-8).  The  resulting  small-signal  state  space  representation  is: 
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(3-11) 
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By  taking  the  determinant  of  the  system  matrix  and  setting  it  equal  to  zero,  a  third-order 
characteristic  polynomial  is  obtained.  The  eigenvalues  of  this  polynomial  determine  the 
closed-loop  pole  positions.  As  a  result,  proper  selection  of  the  control  gains  will  give  the 
desired  pole  locations. 

C.  STATE  DIFFERENCE  IMPLEMENTATION  IN  HARDWARE 


1.  Basic  Implementation 

Implementation  of  Equation  (3-10)  required  generating  the  perturbation  terms, 

v0  (0  and  iL  (t) ,  using  state  difference  methods.  In  this  approach,  v0  (t)  and  iL  (t)  are 

represented  by  difference  terms.  For  example,  the  output  voltage  perturbation  is  given 
by 


vo(0  =  vc-vre/  (3-12) 

where  vre/  is  the  desired  output  voltage  of  the  buck  converter.  As  can  be  seen  from 


Equation  (3-12),  v0(t)  is  zero  in  the  steady  state.  If  a  load  transient  occurs,  however,  vc 
changes  and  produces  a  difference  error  that  is  used  in  Equation  (3-10)  to  affect  the  duty 

cycle  perturbation  term.  Likewise,  the  iL  (/)  difference  term  is 

h(t)  =  i0-iL  (3-13) 
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where  i0  is  output  load  current  and  iL  is  inductor  current.  Substituting  Equation  (3-12) 
and  Equation  (3-13)  into  Equation  (3-10)  results  in  the  state  difference  equation: 

d(t)  =-hi  (i0 -iL)~  hv  (vc  -  vref )  -  h„  J(vc  -  vre/  )d%  (3-14) 

The  development  of  Equation  (3-13)  raises  an  important  point.  Unlike  Equation 
(3-12),  Equation  (3-13)  does  not  employ  the  use  of  a  reference  term  to  act  as  a  baseline 
from  which  to  measure  perturbations.  Rather,  to  generate  difference  errors,  Equation 
(3-13)  relies  on  the  fact  that  inductor  current  will  lag  output  current  during  a  load 
transient.  If  inductor  current  is  an  average  value,  then  Equation  (3-13)  represents  the 
perturbation  in  average  inductor  current  as  developed  in  Equation  (3-10).  Digital 
controller  implementation  allows  the  user  to  meet  this  requirement  as  will  be  shown  in 
Chapter  IV.  Analog  controller  implementation,  however,  does  not  provide  an  adequate 
method  for  obtaining  the  true  average  inductor  current  since  filtering  introduces 
additional  poles  and  unwanted  time  delays.  As  a  result,  Equation  (3-13)  actually 
implements  an  instantaneous  state  difference  for  analog  control.  This  point  raises 
important  gain  selection  issues  outlined  in  Chapter  IV. 

2.  Feedforward  Gain 

Up  to  this  point,  the  small-signal  duty  cycle  algorithm  has  been  determined.  To 
drive  the  switch,  however,  an  expression  for  the  overall  duty  cycle  is  required: 

d(t)  =  D  +  d(t)  (3-15) 
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Equation  (3-15)  states  that  the  overall  duty  cycle  is  made  up  of  a  large-signal  component, 

D,  and  a  small-signal  component,  d(t) ,  which  was  given  in  Equation  (3-14).  An 
expression  must  now  be  developed  for  D. 

Manipulating  Equation  (2-1)  from  Chapter  II,  it  is  apparent  that  D  is  the  ratio  of 
output  voltage  to  input  voltage.  This  ratio  could  be  set  as  a  constant  in  the  controller  and 

methodically  added  to  d(t)  every  switching  cycle,  but  the  full  effectiveness  of  algorithm 
would  not  be  achieved.  In  order  to  enhance  the  control  algorithm,  the  large-signal  duty 
cycle  term  is  implemented  as 

r>  Vre/ 

d"-7Q  P-i«) 

where  vre/is  the  desired  output  voltage  of  the  buck  converter,  and  e(t)  is  the  time- varying 
input  voltage. 

Equation  (3-16)  illustrates  an  important  point.  A  feedforward  gain  is  introduced 
by  using  e(t)  in  the  calculation  of  Equation  (3-16)/  This  gain  compensates  for  the  input 
voltage  perturbations  neglected  in  the  development  of  the  small-signal  duty  cycle 
algorithm.  Since  d(t)  can  instantaneously  adjust  to  input  voltage  changes  due  to  the 
feedforward  gain,  it  acts  to  filter  unwanted  input  noise  from  the  output  of  the  buck 
converter.  The  details  of  this  point  are  expanded  in  Chapter  IV.  Incorporating  Equations 
(3-16)  and  (3-14)  into  Equation  (3-15),  the  total  duty  cycle  algorithm  is  given  by: 

d(t)  =  Dss  -  hi (iQ -iL)~  hv (vc  -  vref )  -  hn  J(vc  -  vref  )d£  (3-17) 
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3. 


Source  Buck  Converter  House  Curve 


With  Equation  (3-17)  in  hand,  a  controller  simulation  and  implementation  can  be 
undertaken  for  an  isolated  buck  converter.  Parallel  operations,  however,  will  be  required 
of  the  Source  Buck  Converters  illustrated  in  Figure  (2-2).  Paralleling  offers  the 
advantages  of  supplying  higher  loads  and  enhancing  reliability  of  the  system.  Before  this 
step  can  be  accomplished,  Equation  (3-17)  must  be  modified  to  allow  load  sharing 
between  buck  converters.  To  accomplish  this  action,  a  house  curve  is  built  into  the 
system  control  and  is  depicted  in  Figure  (3-2).  As  illustrated  in  the  slope  of  Figure  (3-2), 
every  three  amp  increase  in  load  current  causes  the  output  voltage  to  drop  by  one  volt. 


Over  the  no-load  to  full-load  range,  output  voltage  drops  ten  volts.  The  slope  of  the 
house  curve  was  selected  to  provide  an  adequate  load  sharing  profile  while  maintaining 
output  voltage  sufficiently  high  to  supply  system  loads. 
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The  house  curve  is  added  by  introducing  the  load  current,  scaled  by  the  desired 
slope  of  the  house  curve,  into  the  voltage  and  integral  gain  terms  of  Equation  (3-17). 

The  house  curve  could  also  be  introduced  into  the  feedforward  gain  term.  Its  inclusion, 
however,  makes  the  overall  duty  cycle  so  responsive  to  load  current  changes  that 
unacceptable  voltage  and  current  overshoots  occur  during  transient  load  response. 
Equation  (3-18)  is  the  final  form  of  the  analog  control  algorithm: 

<*(')=  D„-h,(i0  -v*  +^j) -h,  j(”c  -V  (3-18) 

This  result  will  be  utilized  in  Chapter  IV  to  simulate  the  analog  control  of  the  Source 
Buck  Converter.  For  use  with  the  Load  Buck  Converter,  the  house  curve  is  removed 
from  Equation  (3-18),  and  it  is  modified  into  the  following  digital  form: 

=  Dss(n)  -  hj (iQ (»)  -  iL  (nj)  -  hvvenor  (n)  -  h„vint  (n)  (3-19) 

where 

T 

Vin»  00  =  {v error  00  +  V error  (»  ~  I))  +  Vta  («  -  1)  (3-20) 

and 

V error  (»)  =  Vc(«)  -  Vfe/  (ft)  (3-21) 

D.  SUMMARY 

This  chapter  detailed  the  development  of  state  difference  control  algorithms  for 
the  Source  and  Load  Buck  Converters.  It  was  shown  that  proper  selection  of  the  control 
gains  will  give  the  desired  pole  locations.  As  a  result,  the  next  step  in  the  system  design 
is  directed  towards  determining  the  gain  magnitudes  utilizing  computer  simulations.  As 
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Chapter  IV  will  show,  this  process  involves  many  tradeoffs  to  achieve  the  required 
system  transient  response. 
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IV.  CONTROLLER  GAIN  SELECTION 


A.  ANALOG  CONTROLLER  GAIN  SELECTION 

1.  Computer  Modeling 

Controller  gains  and  operational  characteristics  were  established  using  computer 
modeling.  First-order  gain  selection  was  done  using  Simulink,  MATLAB’s  dynamic 
system  modeling  tool.  Detailed  analysis  was  then  performed  using  models  developed  in 
the  Advanced  Continuous  Simulation  Language  (ACSL). 

Simulink,  in  combination  with  supporting  MATLAB  m-files,  provided  the  means 
for  iteratively  manipulating  pole  placement  and  fine  tuning  gain  selection.  The  closed- 
loop  Source  Buck  Converter  model  is  pictured  in  Appendix  (B).  The  open-loop  plant  is 
represented  by  the  state  space  averaged  ABCD  matrices.  It  should  be  noted  that  the  state 
space  averaging  process  eliminates  the  switch  by  “averaging”  the  modeling  equations  for 
the  two  topologies.  As  a  result,  characteristics  like  ramping  inductor  current  are  seen  as  a 
dc  average.  The  main  control  algorithm  developed  in  Chapter  III  “closes  the  loop”. 
Utilizing  this  model  and  the  three  MATLAB  m-files  provided  in  Appendix  (B),  Source 
Buck  Converter  gain  selection  performance  was  measured  by  running  a  10%  — »  100%  — » 
10%  load  transient  with  a  hard  input  voltage  source. 

Following  initial  development  in  Simulink,  the  controller  gains  were  evaluated 
under  the  same  10%  — >  100%  — >  10%  load  transient  in  the  detailed  ACSL  model.  Built 
on  the  foundation  of  FORTRAN,  ACSL  allows  modeling  of  the  switched  Source  Buck 
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Converter.  As  a  result,  the  “details”  of  system  performance  are  brought  to  the  surface. 
Model  assumptions  included: 

•  instantaneous  switching 

•  switch  and  diode  losses  modeled  as  voltage  drops 

•  ideal  passive  components 

The  code  for  the  ACSL  simulation  is  provided  in  Appendix  (C). 

2.  Source  Buck  Converter  Transient  Specifications 
Specifications  for  the  Source  Buck  Converter  did  not  include  requirements  for 
system  transient  response.  In  a  general  sense,  transient  response  should  have  short  rise 
and  settling  times  and  little  or  no  overshoot.  In  addition,  the  loop  speed  of  the  controller 
should  be  such  that  low-frequency  (<  360  Hz)  ac  ripple  is  filtered  from  the  input.  This 
general  system  performance  has  been  quantified  in  Table  (4-1).  The  numbers  shown  are 
a  “first  cut”  on  adequate  system  response. 


Figure  of  Merit 

Value 

rise  time  (t,) 

<  1  msec 

2%  settling  time  (ts) 

<  5  msec 

maximum  percent  overshoot  (MJ 

■  okln  A  1  _ HP _ ■ 

5% 

Table  4-1,  System  Transient  Response  Figures  of  Merit 


3.  Closed  Loop  Model 

Before  pole  placement  was  done,  the  Source  Buck  Converter  small-signal  closed- 
loop  block  diagram  was  developed.  Figure  (4-1)  gives  the  final  result.  The  picture  shows 
that  two  effects  can  be  studied.  The  first  involves  analysis  of  the  small-signal  output 

voltage  response  due  to  d .  This  step  is  done  assuming  input  voltage  is  constant. 
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Figure  4-1,  Source  Buck  Converter  Small-Signal  Closed-Loop  System 


As  a  result,  De  is  set  to  zero.  Derivation  of  the  closed-loop  transfer  function  for  this 
model  resulted  in  Equation  (4- 1 ) : 


Vc 


V  ref 


LC 


(< Ch,s 2  +hvs  +  h„ ) 


s 3  + 


1  Eh 


RC 


s2  + 


r\  +  EhA  Eh 

s  +  ■ 


V  LC 


(4-1) 


LC 


The  second  effect  involves  analysis  of  the  input  voltage  ripple  response.  This  step  is 


done  assuming  the  duty  cycle  is  constant.  As  a  result,  vre/  is  set  to  zero.  Derivation  of 
the  closed-loop  transfer  function  for  this  model  resulted  in  Equation  (4-2): 


vc 


D 

LC 


s 3  + 


1  Eh)  2  fl  +  Eh 

VRC+  LJS  H  LC 


s  +  ■ 


E_K 

LC 


(4-2) 


As  shown  below,  both  of  these  cases  will  be  used  to  determine  the  overall  closed-loop 
performance  of  the  Source  Buck  Converter. 


4.  Pole  Placement  And  Gain  Selection 

To  determine  pole  placement,  the  Bessel  Prototype  model  given  in  reference  [6] 
was  employed.  The  poles  for  a  third-order  system  are  given  by 
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s  =  -0.9420wo,  (-0.7455  ±  0.71 12i)w0 


where  w0  is  the  desired  closed-loop  bandwidth.  With  the  poles  established,  a  third-order 
polynomial  can  be  formed,  and  term-by-term  matching  can  be  done  with  the  denominator 
of  Equation  (4-1)  to  find  the  controller  gains. 

Selection  of  w0,  and  thus  the  gains,  was  set  by  three  factors.  First,  w0  must  be 
sized  to  ensure  the  pole  farthest  in  the  left  half  plane  is  at  least  a  factor  often  smaller  than 
the  radian  switching  frequency  to  prevent  unwanted  controller  interactions.  Second, 
selection  of  w0  should  not  require  an  excessive  duty  cycle  control  effort.  This  effect 
causes  increased  noise  in  the  voltage  output  due  to  a  beating  action  set  up  between  the 
faster  duty  cycle  response  and  the  slower  output  capacitor  voltage  response.  Third,  w0 
must  be  large  enough  to  produce  the  required  response  times  of  Table  (4-1). 


Based  on  these  criteria,  the  final  pole  locations  and  controller  gains  were 
determined  using  w0  =  3250  rad/sec.  The  results  are  documented  in  Table  (4-2). 


Closed-Loop  Pole  Locations 

s  =-3062, -2423  ±23 lli  1 

Proportional  Voltage  Gain 

o 

o 

II 

Integral  Voltage  Gain 

h„  =  26.09 

Proportional  Current  Gain 

'T'  l_l  _  A  C\ _  Tk  i  ^ 

hj  =  0.015  ] 

Table  4-2,  Source  Buck  Converter  Closed-Loop  Poles  and  Gains 


Substituting  these  gain  values  into  Equation  (4-1),  the  system  closed-loop  frequency 
response  shown  in  Figure  (4-2)  was  produced.  From  the  magnitude  response  it  can  be 
seen  that  frequencies  beyond  7000  rad/sec  (~1 1 14  Hz)  are  attenuated  by  the  20  dB/decade 
rolloff.  Up  to  this  point,  the  response  shows  little  or  no  attenuation.  In  fact,  from  2000 
rad/sec  to  7000  rad/sec,  a  maximum  gain  of  2  dB  occurs.  The  input  ripple  response  of 
Figure  (4-3),  however,  compensates  for  this  lack  of  attenuation.  This  plot  shows  a 
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A  j  A 

Figure  4-2,  Source  Buck  Converter  Closed-Loop  Frequency  Response:  vcf  vre/ 


Figure  4-3,  Source  Buck  Converter  Closed-Loop  Frequency  Response:  vc/<? 
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minimum  attenuation  of  20  dB  in  the  360  Hz  frequency  range.  As  a  result,  the  combined 
effect  of  both  responses  yields  proper  Source  Buck  Converter  closed-loop  operations. 

The  MATLAB  transient  response  is  given  in  Figure  (4-4).  Analysis  of  Figure  (4-4)  and 
the  MATLAB  code  results  shows  that  all  the  requirements  of  Table  (4-1)  were  met. 

Next,  the  transient  response  was  verified  using  the  detailed  ACSL  simulation. 


Substituting  the  Table  (4-2)  gains  into  the  ACSL  model  produced  the  transient 


response  pictured  in  Figure  (4-5).  The  most  notable  feature  of  the  ACSL  graphs  is  the 
influence  of  detailed  switching  on  inductor  current  and  duty  cycle.  The  inductor  current 
response  was  expected.  The  duty  cycle  response,  on  the  other  hand,  was  not.  In  fact, 
results  close  to  that  of  the  duty  cycle  plot  in  Figure  (4-4)  were  anticipated.  The  reason  for 
the  duty  cycle  band  in  Figure  (4-5)  goes  back  to  the  fact  that  the  inductor  current 
perturbation  term  in  Equation  (3-18)  is  not  a  true  average.  Rather,  instantaneous 
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Figure  4-5,  ACSL  Source  Buck  Converter  Closed-Loop  Transient  Response 
inductor  current  and  load  current  are  subtracted  to  produce  the  perturbation  term.  In 
steady  state,  this  difference  is  the  inductor  ac  current  ripple.  As  a  result,  this  ripple  is 
introduced  into  the  duty  cycle  calculation,  and  the  duty  cycle  waveform  of  Figure  (4-5)  is 
produced.  The  end  result  of  this  effect  is  the  introduction  of  low  amplitude  (<  0.1  V) 

20  kHz  harmonics  into  the  output.  Therefore,  it  is  imperative  to  minimize  the  magnitude 
of  h;  while  still  meeting  the  response  of  Table  (4-1).  Bench  testing  showed  that  a  value  of 
1^  less  than  0.02  gave  acceptable  results. 

5.  Effect  of  Output  Capacitor  Sizing 

To  this  point,  pole  placement  and  gain  selection  have  been  done  utilizing  an 
output  inductance  of  760  pH  and  an  output  capacitance  of  400  pF.  It  has  been  shown, 
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however,  that  output  capacitance  can  be  as  low  as  10  pF  while  still  satisfying  the  1% 
output  ripple  requirement.  As  a  result,  justification  for  the  400  pF  choice  is  required. 

Using  the  Bessel  Prototype  model  described  above,  system  transient  response  was 
analyzed  at  additional  output  capacitance  values  of  10  pF,  140  pF,  230  pF,  and  400  pF. 
A  value  of  w0=  3250  rad/sec  was  used  for  all  cases  to  maintain  the  current  gain  less  than 
or  equal  to  its  present  value  based  on  the  discussion  in  the  previous  paragraph. 

Table  (4-3)  gives  the  gains  for  each  of  the  output  capacitances  mentioned  above.  Figure 
(4-6)  summarizes  the  transient  responses. 


Parameter 

C  =  10  pF 

C  =  140  pF 

C  =  230  pF 

C  =  400  pF 

hi 

0.0131 

0.015 

r  0.015 

0.015 

hv 

-0.0021 

0.0043 

0.0087 

0.017 

K 

0.65 

9.13 

15.00 

26.09 

1  able  4-3,  Capacitor  Sizing  Effects  on  the  Source  Buck  Converl 

ter  Gains 

Vc  Response:  C  =  10  uF 


time  (msec) 


Figure  4-6,  Source  Buck  Converter  Transient  Response  as  a  Function  of  Capacitor  Size 
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From  the  data  provided,  it  can  be  seen  that  the  10  jiF  case  was  unacceptable.  All  of  the 
remaining  cases,  however,  were  well  within  the  Table  (4-1)  requirements.  Based  on 
j  overshoot  and  settling  characteristics,  however,  the  400  pF  capacitor  was  used  in  the 

hardware  implementation.  This  selection  was  made  to  prevent  undesirable  buck 
converter  interactions  during  parallel  operations  and  to  ensure  adequate  system  response 
at  load  transient  testing  above  the  10%  ->  100%  10%  scenario  used  in  this  analysis. 

B.  DIGITAL  CONTROLLER  GAIN  SELECTION 

1.  Computer  Modeling 

As  with  the  analog  controller,  the  digital  controller  was  modeled  using  Simulink 
and  ACSL.  The  simulation  limitations  discussed  above  also  apply  to  this  section.  In  the 
digital  simulation,  the  continuous-time  controller  inputs  were  sampled  once  each 
switching  period  using  a  zero-order  hold  and  applied  to  the  digital  control  algorithm 
given  by  Equation  (3-19).  The  calculated  duty  cycle  is  then  applied  to  the  IGBT  gate 
drive  on  the  next  switching  period.  The  Simulink  block  diagram  and  MATLAB  m-files 
are  contained  in  Appendix  (B).  The  code  for  the  ACSL  simulation  is  provided  in 
Appendix  (C). 

2.  Pole  Placement  and  Gain  Selection 

The  design  parameters  of  Table  (4-1)  and  the  simulation  sequence  employed 
above  were  also  used  for  the  digital  controller  design.  In  addition,  the  limitations  on  w0 
also  hold  here.  Unlike  the  analog  case,  however,  current  gain  selection  does  not  play  as 
significant  a  role.  This  fact  was  due  to  the  use  of  a  true  average  for  the  inductor  current 
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perturbation  term  in  Equation  (3-19).  Since  the  analog  inputs  to  the  digital  controller 
are  sampled  once  each  switching  period,  the  inductor  current  sample  becomes  an  average 
value  as  long  as  the  sample  point  occurs  at  the  same  place  during  each  switching  period. 
Of  course  the  ramping  nature  of  inductor  current  may  cause  the  sampled  average  to  be 
slightly  different  from  the  actual  average  depending  where  on  the  ramp  the  sample  is 
taken.  Any  offset,  however,  is  compensated  for  by  the  digital  integrator.  This  argument 
is  equally  applicable  to  the  case  where  multiple  samples  are  taken  and  averaged  during 
the  switching  period. 

3.  Summary  of  Digital  Controller  Gain  Selections 

Table  (4-4)  provides  a  summary  of  final  gain  selections  for  the  Load  Buck 
Converter  controller.  The  design  centered  on  using  a  400  pF  capacitor  on  the  output  with 
w0  =  1900  rad/sec.  The  value  of  w0  was  set  based  on  a  maximum  duty  cycle  limitation  of 
0.95  for  the  Universal  Controller  detailed  in  reference  [7].  Since  the  controller  was  built 
for  dual  Auxiliary  Resonant  Commutated  Pole  (ARCP)  Inverter  and  Buck  Converter 
operations,  the  duty  cycle  must  incorporate  dead  time  to  prevent  a  direct  short  between 
ARCP  switches  during  the  switching  interval.  As  a  result,  duty  cycle  is  limited  to  0.95. 
The  w0  selection  prevents  the  duty  cycle  from  exceeding  this  limit  so  that  linear  operation 
is  maintained.  Overall,  Table  (4-1)  transient  response  is  achieved.  See  the  analog 


Closed-Loop  Pole  Locations 

s  =-942,  -746  ±711i 

Proportional  Voltage  Gain 
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Table  4-4,  Load  Buck  Converter  Closed-Loop  Poles  and  Gains 
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controller  case  above  for  the  detailed  gain  and  capacitor  selection  procedure.  Figure  (4-7) 


documents  the  ACSL  transient  load  response. 


Figure  4-7,  ACSL  Load  Buck  Converter  Closed-Loop  Transient  Response 


C.  SUMMARY 

With  the  controller  gains  established  for  the  Source  and  Load  Buck  Converters, 
simulation  of  the  PEBB  Network  can  begin.  This  analysis  will  be  used  to  identify 
stability  criteria  for  the  network.  With  this  information  in  place,  interesting  PEBB 
Network  component  interactions  can  be  investigated. 
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V. 


SIMULATION  AND  HARDWARE  RESULTS 


A.  BACKGROUND 

With  the  proper  analog  and  digital  controller  gains  established,  detailed  PEBB 
Network  simulations  were  developed  and  conducted  in  ACSL.  Representative  models, 
identified  as  “modes”  in  this  chapter,  were  established  based  on  seven  basic  PEBB 
Network  hardware  configurations.  Model  assumptions  included: 

•  rectified  3-phase  power  input  to  Source  Buck  Converters 

•  instantaneous  switching 

•  no  switching  noise 

•  switch  and  diode  losses  modeled  as  voltage  drops 

•  inductor  and  capacitor  equivalent  series  resistance  (ESR)  modeled 

•  DC  bus  resistive-inductive  link  (RL-link)  effects  not  modeled 

The  RL-link  effect  was  not  considered  due  to  a  tradeoff  between  simulation  run  time  and 
the  simulation  time  step.  The  time  step  could  be  made  small  enough  to  model  the  RL- 
links,  but  the  simulation  run  time  would  have  been  unnecessarily  large.  In  addition,  the 
RL-link  contribution  to  the  overall  system  stability  is  negligible  for  the  components 
considered  based  upon  the  small  size  of  the  line  inductance  compared  to  the  much  larger 
inductance  value  of  the  Load  Buck  Converter  input  filter. 

For  each  of  the  seven  modes  described  in  this  chapter,  two  transient  analysis  cases 
were  considered.  The  first  case  investigated  a  step  change  in  load  resistance  on  the 
output  of  the  Source  Buck  Converter.  Load  Buck  Converter  output  response  was  then 
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observed  to  identify  if  its  controller  could  regulate  the  desired  output  voltage  and  current. 
No  noticeable  effects  or  stability  issues  were  documented  for  any  of  the  seven  topologies 
(modes),  and  therefore,  this  transient  study  was  not  considered  further.  The  second 
transient  case  performed  a  load  change  on  the  output  of  the  Load  Buck  Converter  and/or  a 
conventional  inverter.  Overall  system  response  was  then  observed.  The  results  for  this 
case  are  documented  below  for  each  of  the  seven  modes.  Representative  ACSL  code  is 
provided  for  Mode  7  operation  in  Appendix  (C). 

B.  MODE  1  OPERATION 


Mode  1  is  defined  as  a  Source  Buck  Converter  operating  in  series  with  a  Load 
Buck  Converter.  Figure  (5-1)  illustrates  the  overall  system  model.  The  only  difference 


component-wise  between  Figure  (5-1)  and  the  ACSL  code  resulted  from  the  need  to 
represent  the  input  filter  to  the  Source  Buck  Converter  differently  in  the  simulation.  Due 
to  an  algebraic  loop  formed  between  the  ACSL  main  program,  the  buck  converter  macro 
and  the  analog  controller  macro,  the  input  filter  capacitor  and  its  associated  ESR  had  to 
be  modeled  by  a  capacitor  in  parallel  with  an  equivalent  resistor.  As  a  result,  the 

Rcesrjnl’  component  in  Figure  (5-1)  is  represented  by  a  2.47  kQ  parallel  resistance  in 
the  ACSL  code. 
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System  transient  analysis  focused  on  investigating  a  step  load  change  from  300  W 
to  3  kW  on  the  output  of  the  Load  Buck  Converter.  Figure  (5-2)  documents  the  results. 

As  can  be  seen,  the  output  voltage  and  current  of  the  Source  Buck  Converter  (Vout  and 
ioutl)  have  a  somewhat  more  pronounced  oscillation  at  the  transient  edges  before  settling 
to  steady-state  values.  This  result  is  due  to  a  “ringing”  effect  set  up  between  the  output  of 
the  Source  Buck  Converter  and  the  input  LC  filter  of  the  Load  Buck  Converter.  In  fact, 
this  effect  is  clearly  discernible  in  ioutl  in  the  steady  state.  The  transient  performance  of 
the  Load  Buck  Converter  is  the  same  as  that  predicted  by  Figure  (4-7)  in  Chapter  IV. 

With  system  stability  demonstrated,  a  hardware  validation  study  was  conducted. 
Figure  (5-3)  shows  the  actual  system  response  for  a  300  W  to  2.5  kW  load  change.  The 
reason  for  the  difference  between  the  simulated  3  kW  maximum  load  value  and  the  actual 
2.5  kW  load  value  was  based  upon  the  sizing  of  test  load  banks  available  in  the 
laboratory.  Comparison  of  the  simulation  and  hardware  trials  demonstrates  reasonably 
good  correlation  between  predicted  and  actual  system  response.  Note  that  the  voltages  in 
Figure  (5-3)  are  ac  coupled.  The  major  difference  between  the  plots  comes  from  the 
switching  harmonics  and  noise  of  the  actual  system  operation.  The  principle  source  of 
the  noise  comes  from  the  operation  of  the  IGBT  switches. 

C.  MODE  2  OPERATION 

Mode  2  is  defined  as  a  Source  Buck  Converter  operating  in  series  with  an 
inverter.  Figure  (5-4)  illustrates  the  overall  system  model.  A  hard-switched  inverter 
was  used  to  model  the  ARCP  unit  so  that  simulation  times  could  be  kept  to  a  reasonable 
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length.  Since  the  gross  system  dynamics  between  an  ARCP  and  a  conventional  inverter 
are  nearly  the  same  this  assumption  was  reasonable  to  make.  The  inverter  was  controlled 
by  a  current  feedback  controller  implemented  in  the  stationary  reference  frame.  A 
complete  discussion  on  the  derivation  of  this  model  is  given  in  reference  [8].  As  with 
Mode  1  operation,  the  input  filter  capacitor  and  its  associated  ESR  had  to  be  modeled  as 
an  equivalent  parallel  combination  in  the  ACSL  code. 


System  transient  analysis  focused  on  investigating  a  step  change  in  load  from 
2.5  kW  to  6  kW  at  the  inverter  output.  This  was  performed  by  changing  the  amplitude 
values  of  the  reference  phase  currents  from  10  A  to  24  A.  Figure  (5-5)  documents  the 
results.  As  can  be  seen,  significant  switching  harmonics  have  been  introduced  into  the 
output  voltage  and  current  responses  of  the  Source  Buck  Converter.  In  addition,  this 
effect  is  aggravated  at  higher  levels  of  inverter  loading.  The  output  of  the  inverter  is  as 
expected  for  a  current  controlled  unit. 


40 


One  point  to  note  is  that  the  inverter  current  response  does  not  achieve  the  desired 
steady-state  values.  This  result  stems  from  the  use  of  stationary  reference  frame  control. 
Although  it  is  easy  to  implement,  reference  [8]  shows  that  steady-state  errors  of  11%  can 
be  expected.  A  hardware  validation  study  was  not  conducted  for  this  case  because  the 
ARCP  closed-loop  control  has  not  been  fully  developed  at  this  time.  It  is  anticipated  that 
future  research  efforts  will  provide  the  hardware  study  results  to  validate  these 
waveforms. 


D.  MODE  3  OPERATION 

Mode  3  is  defined  as  a  Source  Buck  Converter  supplying  both  a  Load  Buck 
Converter  and  an  inverter.  Figure  (5-6)  illustrates  the  overall  system  model.  The  same 
inverter  model  described  in  Mode  2  operation  was  used  for  this  case.  As  with  Mode  1 
operation,  the  input  filter  capacitor  and  its  associated  ESR  had  to  be  modeled  as  an 
equivalent  parallel  combination  in  the  ACSL  code. 
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Transient  analysis  focused  on  investigating  system  response  to  the  load  changes 
discussed  in  Mode  1  and  Mode  2  operation.  The  step  change  for  the  Load  Buck 
Converter  occurred  at  45  msec  while  the  step  change  in  inverter  loading  was  specified  at 
125  msec.  Figure  (5-7)  documents  the  results.  The  plots  show  a  system  response  that  is 
the  composite  of  the  results  achieved  in  the  previous  two  modes  of  operation.  As  with 
Mode  2  operation,  a  hardware  validation  study  was  not  conducted  since  the  ARCP 
closed-loop  control  algorithm  is  still  under  development. 

E.  MODE  4  OPERATION 

Mode  4  is  defined  as  parallel  Source  Buck  Converter  operations.  Figure  (5-8) 
illustrates  the  overall  system  model.  As  described  in  Mode  1  operation,  the  input  filter 
capacitors  and  their  associated  ESRs  had  to  be  modeled  as  equivalent  parallel 
combinations  in  the  ACSL  code. 

Transient  analysis  for  this  case  focused  on  investigating  a  load  change  from 
2.25  kW  to  18  kW  and  observing  system  response.  Figure  (5-9)  documents  the  results. 
From  the  plots,  it  can  be  seen  that  the  parallel  units  share  the  total  load,  iout,  equally.  In 
addition,  the  output  voltage  response  covers  the  full  range  of  house  curve  operation  as 
expected. 

An  interesting  occurrence  to  note  in  the  individual  buck  converter  current 
responses,  ioutl  and  iout2,  is  the  “mirror”  effect  seen  at  the  transient  edges.  This  effect  is 
best  observed  in  the  system  transient  at  the  85  msec  mark.  At  this  point,  ioutl  settles  to 
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Figure  5-8,  PEBB  Testbed  Mode  4  Operation 


its  final  steady-state  value  following  a  slight  undershoot.  The  response  for  iout2,  on  the 
other  hand,  settles  to  its  final  steady-state  value  after  experiencing  an  overdamped  effect. 
When  observed  on  an  expanded  scale,  the  currents  are  mirror  responses  of  each  other. 

This  result  occurs  due  to  the  difference  in  the  output  inductor  sizing  of  the  Source  Buck 
Converter  units.  In  essence,  a  small  AC  current  is  established  and  reflects  between  the 
paralleled  capacitors  based  upon  the  difference  in  the  inductor  current  waveforms. 

With  system  stability  demonstrated  and  the  transient  excursion  of  the  variables 
well  within  device  limitations,  a  hardware  validation  study  was  conducted  for  Mode  4 
operation.  Figure  (5-10)  shows  the  actual  system  response  for  a  2.25  kW  to  7.5  kW  load 
change.  Although  the  overall  transient  was  smaller  in  magnitude,  a  general  comparison 
of  the  simulation  and  hardware  trials  demonstrates  reasonably  good  correlation  between 
predicted  and  actual  system  response.  As  with  the  Mode  1  hardware  validation,  the  major 
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difference  between  the  simulated  and  actual  results  is  due  to  the  presence  of  switching 
harmonics  and  EMI  noise  in  actual  system  operation. 

The  current  “mirror”  effect  discussed  previously  is  demonstrated  on  an  expanded 
scale  in  Figure  (5-11)  for  actual  circuit  operation.  It  can  be  seen  that  the  effect  occurs 
during  steady  state  as  well  as  transient  operations.  The  steady-state  effect  was  not  as 
noticeable  in  the  simulation  because  switching  noise  was  not  modeled.  Follow-on 
research  may  consider  incorporating  a  Gaussian  white  noise  source  into  the  sensed 
variables  in  the  simulation. 

F.  MODE  5  OPERATION 

Mode  5  is  defined  as  two  parallel  Source  Buck  Converters  operating  to  supply  one 
Load  Buck  Converter.  Figure  (5-12)  illustrates  the  overall  system  model.  As  described 
in  Mode  1  operation,  the  input  filter  capacitors  and  their  associated  ESRs  had  to  be 
modeled  as  equivalent  parallel  combinations  in  the  ACSL  code. 

System  transient  analysis  focused  on  investigating  a  step  change  in  resistive  load 
from  300  W  to  3  kW  on  the  output  of  the  Load  Buck  Converter.  Figure  (5-13) 
documents  the  simulation  results.  As  can  be  seen,  the  output  voltage  and  current  of  the 
parallel  Source  Buck  Converters  have  a  somewhat  more  pronounced  oscillation  at  the 
transient  edges  before  settling  to  their  steady-state  values.  This  result  is  due  to  a 
“ringing”  effect  set  up  between  the  output  of  the  Source  Buck  Converter  and  the  input  LC 
filter  of  the  Load  Buck  Converter.  During  load  transients,  this  effect  dominates  over  the 
current  “mirror”  effect  documented  in  Mode  4  operation.  This  result  is  shown  in 


47 


i°ut2  (amps)  ioutl  (amps) 


Figure  (5-13)  by  observing  that  the  transient  portions  of  the  output  currents,  ioutl  and 
iout2,  are  “in  phase”  with  one  another  during  load  changes. 

With  system  stability  demonstrated,  a  hardware  validation  study  was  conducted 
for  Mode  5  operation.  Figure  (5-14)  shows  the  actual  system  response.  From  the  plots,  it 
can  be  seen  that  the  actual  Load  Buck  Converter  output  correlated  with  the  simulation 
results.  Actual  current  values  out  of  the  parallel  Source  Buck  Converters  were  about 
2.5  A  higher  than  shown  in  the  simulation  results  because  a  25  Q  static  load  was  used 
during  hardware  validation.  In  the  ACSL  simulation,  this  value  was  40  Q.  The  voltage 
output  response  of  the  parallel  Source  Buck  Converter  is  harder  to  compare  due  to  the 
switching  noise  present  on  the  waveform.  Note  that  all  voltage  responses  are  ac  coupled. 

The  “ringing”  effect  versus  the  current  “mirror”  effect  is  shown  in  Figure  (5-15). 
This  plot  shows  how  the  mirror  effect  dominates  dining  steady-state  operation  but  shifts 
“in  phase”  during  load  changes.  Overall,  comparison  of  the  simulation  and  hardware 
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Figure  5-14,  Mode  5  Hardware  Results:  300  W  2.5  kW  Load  Buck  Transient 
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trials  demonstrates  reasonable  correlation  between  predicted  and  actual  results.  As  with 
the  previous  hardware  studies,  switching  noise  constitutes  the  major  difference  between 
simulated  and  actual  operations. 


G.  MODE  6  OPERATION 

Mode  6  is  defined  as  two  parallel  Source  Buck  Converters  supplying  an  inverter. 
Figure  (5-16)  illustrates  the  overall  system  model.  The  same  inverter  model  described  in 
Mode  2  operation  was  used  for  this  case.  As  described  in  Mode  1  operation,  the  input 
filter  capacitors  and  their  associated  ESRs  had  to  be  modeled  as  equivalent  parallel 
combinations  in  the  ACSL  code. 
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System  transient  analysis  focused  on  investigating  a  step  change  in  inverter  load 
from  2.5  kW  (10  A  peak  current  per  phase)  to  6  kW  (24  A  peak  current  per  phase)  and 
observing  system  response.  Figure  (5-17)  documents  the  results.  As  can  be  seen, 
significant  switching  harmonics  have  been  introduced  into  the  output  voltage  and  current 
responses  of  the  parallel  Source  Buck  Converters.  In  addition,  the  level  of  the  harmonics 
is  increased  with  increased  in  inverter  load.  The  output  of  the  inverter  is  as  expected  for  a 
current  controlled  unit.  As  described  earlier,  the  inverter  current  response  does  not 
achieve  the  desired  steady  state- values  due  to  the  use  of  stationary  reference  frame 
control.  A  hardware  validation  study  was  not  conducted  for  this  case  because  the  ARCP 
closed-loop  control  has  not  been  fully  developed. 


! 
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H.  MODE  7  OPERATION 

Mode  7  is  defined  as  two  parallel  Source  Buck  Converters  supplying  both  a  Load 
Buck  Converter  and  an  inverter.  Figure  (5-18)  illustrates  the  overall  system  model.  The 
same  inverter  model  described  in  Mode  2  operation  was  used  for  this  case.  As 


described  in  Mode  1  operation,  the  input  filter  capacitors  and  their  associated  ESRs  had 
to  be  modeled  as  equivalent  parallel  combinations  in  the  ACSL  code. 


System  transient  analysis  focused  on  investigating  the  transients  discussed  in 
Mode  5  and  Mode  6  operation  and  observing  system  response.  Figures  (5-19)  and  (5-20) 
document  the  results.  The  step  change  for  the  Load  Buck  Converter  occurred  at  45  msec 
while  the  step  change  in  inverter  loading  was  specified  at  125  msec.  These  plots  show  a 
system  response  that  is  the  composite  of  Mode  5  and  Mode  6  results. 
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I.  SUMMARY 

This  chapter  has  included  a  discussion  of  the  ACSL  simulations  performed  for 
seven  basic  PEBB  Network  configurations.  The  results  of  these  simulations  did  not  yield 
any  gross  system  instabilities.  At  most,  only  "interesting”  system  interactions  were 
identified.  In  addition,  hardware  validation  studies  were  conducted  for  Modes  1, 4  and  5. 
The  results  of  these  tests  verified  that  the  ACSL  simulations  correlated  well  with 
hardware  implementation.  The  hardware  validation  results  also  showed  that  switching 
noise  was  a  major  issue  in  actual  circuit  operation.  Chapter  VII  addresses  methods  for 
eliminating  or  minimizing  the  switching  noise  and  suggests  follow-on  ACSL  simulations 
and  hardware  validation  studies. 
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VL  SCHEMATICS 


A.  CIRCUIT  BOARD  BACKGROUND 

1.  Easytrax  Overview 

The  development  of  printed  circuit  boards  (PCBs)  was  done  using  Protel 
Easytrax,  version  2.06.  This  software  allowed  PCB  definition,  layout,  and  component 
interconnection  on  two  layers.  All  PCB  layouts  were  done  so  that  there  were  no  isolated 
ground  sections.  Once  a  PCB  was  developed,  a  netlist  was  compiled.  The  basic  function 
of  the  netlist  is  to  associate  integrated  circuits  with  resistors,  capacitors,  and  diodes  so 
that  connections  can  be  verified.  With  the  PCB  layout  completed,  a  Gerber  file  was  made 
to  direct  the  required  circuit  milling.  The  milling  work  was  done  at  the  Naval 
Postgraduate  School  by  DPI  Daly.  Finally,  hardcopies  of  circuit  traces  and  component 
layouts  were  made  for  future  reference. 

2.  Circuit  Labeling 

For  the  purpose  of  more  easily  locating  components  and  understanding  their 
function  in  the  circuit,  the  following  format  was  chosen: 

The  first  digit  represents  the  type  of  component. 

•  R:  resistor 

•  C:  capacitor 

•  D:  diode 

The  second  digit  represents  the  IC  with  which  the  component  is  best  associated. 
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•  0:  U10 


•  1:  U1 

•  2:  U2,etc. 

The  third  and  fourth  digits  are  used  for  IC  pin  location. 

•  03:  pin  3 

•  18:  pin  18,  etc. 

When  a  prefix  is  used,  it  indicates  an  interface  between  circuits  in  the  overall  system. 
The  prefix  may  have  an  actual  number  associated  with  it  or  just  the  sign. 

•  SI:  Sensor  board  interface  for  buck  converter  1 

•  D#:  Driver  board  interface  for  buck  converter  being  referenced 

•  J25:  25  pin  D-sub  connector  interfacing  controller  with  other  boards 
If  a  suffix  is  used,  then  there  are  generally  more  than  one  of  the  same  components 
connected  to  the  same  pin. 

•  Multiple  components  of  same  type  are  connected  to  same  pin 
where  A  =  1  component,  B  =  2  components,  etc. 

•  G:  Ground 

•  H:  Vcc,  Vdd:  positive  source  voltage 

•  L:  Vss:  negative  source  voltage 

•  Q:  Transistor 

Specific  examples  of  circuit  labeling  include: 

•  C#-J25-9.  Control  board  interface  at  J25  pin  9  for  buck  converter  being 
referenced. 
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•  R1 13C:  Resistor  connected  to  U1  pin  13  and  there  are  at  least  two  more 
resistors  connected  to  that  same  node. 

•  C006:  Capacitor  connected  to  U10  pin  6. 

•  RJ4036Q:  Resistor  connected  to  J40  pin  36  and  a  transistor. 

B.  CIRCUIT  DIAGRAMS  AND  DESCRIPTIONS 

This  section  documents  the  hardware  implementation  of  the  Source  and  Load 
Buck  Converters  as  well  as  their  associated  sensor,  controller,  and  driver  boards.  Each 
subsection  is  layed  out  so  that  it  contains  a  detailed  line  diagram  and  description  of  the 
circuit.  If  further  information  is  needed,  Appendix  (D)  contains  all  the  Easytrax  PCB 
schematics  and  netlists. 

1.  Source  and  Load  Buck  Converter  Topologies 

Figures  (6-1)  and  (6-2)  document  the  layout  of  the  Source  and  Load  Buck 
Converters.  As  seen  from  the  diagrams,  the  topologies  are  identical.  Input/output 
characteristics,  inductor  sizing,  and  the  voltage  sensing  circuits  represent  the  only 
differences.  The  reason  for  the  different  sensing  circuits  goes  back  to  the  choice  of 
analog  control  boards  for  the  Source  Buck  Converters  and  digital  control  boards  for  the 
Load  Buck  Converters.  As  the  diagrams  indicate,  the  analog  controller  requires  a  voltage 
reduction  of  (1/50)  whereas  the  digital  controller  requires  a  voltage  reduction  of  (1/100). 
The  general  flow  of  overall  system  operation  is  as  follows: 

•  Sensed  voltage  from  the  voltage  divider  network  is  sent  to  the  sensor  board 
via  twisted  pair. 
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•  Sensed  currents  are  sent  to  the  sensor  board  via  the  Hall  Effect  sensors. 

•  The  sensor  board,  acting  as  an  isolator,  provides  the  appropriately  scaled 
voltages  and  currents  to  the  analog/digital  control  board. 

•  The  control  board  processes  the  voltages  and  currents  according  to  the 
developed  control  algorithm. 

•  The  generated  driver  signal  from  the  control  board  is  sent  to  the  IGBT  driver 
board  via  twisted  pair  (analog  control)  or  optical  signal  (digital  control). 
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•  The  driver  signal  controls  the  IGBT  switching  frequency  and  duty  cycle. 

With  this  overview  in  hand,  a  more  detailed  description  is  presented  next. 

2.  Sensor  Board 

Pictured  in  Figure  (6-3),  the  sensor  board  provides  isolation  between  the  buck 
converter  and  the  analog/digital  control  board.  Sensed  voltages  from  the  voltage 
divider  network  on  the  buck  converter  are  sent  via  twisted  pair  to  the  input  of  wideband 
AD215  isolation  units.  Configured  as  a  unity  gain  buffers,  the  AD215s  output  the  sensed 
voltages  and  send  them  to  the  control  board  via  twisted  pair.  Current  sensing  is  done  by 
CL50  Hall  Effect  sensors  mounted  directly  on  the  sensor  board.  Instantaneous  currents 
proportional  to  (i^ed/1000)  are  output  by  the  Hall  Effects.  These  currents  are  converted 
to  a  proportional  voltage  by  100  Q  precision  resistors  to  yield  (isense</10).  These  voltages 
are  the  sent  to  the  control  board  via  twisted  pair. 

3.  Analog  Control  Board 

The  implementation  of  the  analog  control  board  can  be  broken  down  into  five 
subsections: 

•  J25  -  25  Pin  Connector  and  Buffer  Stage 

•  Main  Control  Stage 

•  Pulse  Width  Modulation  Stage 

•  Protection  and  Startup  Circuitry  Stage 

•  Power  Supply  Stage 
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A  diagram  and  discussion  of  each  of  the  subsections  is  provided  below.  Appendix  (D) 
contains  detailed  PCB  information  on  the  analog  control  board. 
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Figure  6-3,  Sensor  Board 
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Buffer  Stage 


a. 

From  the  twisted  pair  output  of  the  sensor  board,  the  sensed  and  scaled 
voltages  and  currents  go  to  the  25  pin  D-sub  connector  of  the  control  board  as  shown  in 
Figure  (6-4).  From  J25,  the  inputs  are  buffered  utilizing  a  LM324  quad  op  amp  (labeled 
U2).  Each  buffer  stage  input  contains  an  RC  lowpass  filter  designed  to  remove  high- 
frequency  components  from  the  measured  voltages  and  currents.  The  cutoff  frequency 
for  each  filter  is  at  8  kHz.  The  LM324  outputs  are  fed  to  the  Main  Control  Stage. 

b.  Main  Control  Stage 

The  Main  Control  Stage  is  where  Equation  (2-18)  is  formulated.  The 
output  of  this  stage  is  the  analog  duty  cycle  scaled  by  a  factor  of  ten.  In  addition  to  the 
Buffer  Stage  outputs,  the  signal  representing  the  desired  output  voltage  is  fed  into  the 
Main  Control  from  the  Protection  and  Startup  Circuitry  Stage.  As  shown  in  Figure  (6-5), 
an  AD534  configured  as  a  divider  provides  the  base  duty  cycle,  Dss,  by  dividing  the 
reference  voltage  by  the  input  voltage.  Note  that  this  arrangement  implements  the 
feedforward  term  given  in  Equations  (3-18)  and  (3-19)  of  Chapter  III. 

The  small-signal  control  algorithm  is  performed  using  a  LF347  quad  op 
amp.  The  first  three  op  amps  in  the  package  form  the  proportional  voltage  (Ul,),  integral 
voltage  (Ul2)and  current  response  (Ul3)  terms.  Note  that  a  limited  integrator  has  been 
used  to  generate  the  integral  voltage  term.  This  implementation  prevents  integrator 
windup.  Without  the  limiting  zener  diodes  in  place,  Cl 07  is  free  to  charge  up  during 
extended  transients  to  values  well  outside  the  range  in  which  its  voltage  affects  duty 
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cycle.  As  a  result,  the  voltage  on  Cl 07  locks  out  the  remaining  terms  in  the  duty  cycle 
control  algorithm.  This  lockout  persists  until  Cl 07  come  backs  within  its  operating 
range.  The  final  op  amp  (Ul4)  is  used  to  scale  and  sum  the  outputs  of  the  first  three  op 
amps.  Scaling  is  performed  such  that  the  analog  gain  values  developed  in  Chapter  IV 
were  obtained,  scaled  by  a  factor  of  ten.  Output  feeds  the  Pulse  Width  Modulation  Stage. 


Figure  6-4,  Analog  Control  Board:  J25  -  25  Pin  Connector  and  Buffer  Stage 
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Figure  6-5,  Analog  Control  Board:  Main  Control  Stage 
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c.  Pulse  Width  Modulation  Stage 

Pictured  below  in  Figure  (6-6),  the  Pulse  Width  Modulation  (PWM)  Stage 
is  comprised  of  the  UC3637  chip.  It  receives  the  analog  duty  cycle  from  the  Main 
Control  Stage  and  produces  a  square  wave  driver  signal  with  a  frequency  of  23  kHz  based 
upon  R3 1 8  and  C302.  This  frequency  can  be  calculated  utilizing  formulas  given  in  the 
UC3637  data  sheet  in  Appendix  (A).  Actual  frequency  runs  about  20.5  kHz  due  to 
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interactions  with  the  phase-locked  loop  circuit  described  below.  The  duty  cycle  of  the 
pulses  is  equal  to  the  Main  Control  Stage  analog  duty  cycle  signal  divided  by  ten.  From 

the  Pulse  Width  Modulation  Stage,  the  driver  signal  is  fed  to  a  driver  board  mounted  on 
the  IGBT. 


d.  Phase  Locked  Loop 

Analysis  of  Figure  (6-6)  shows  MASTER  and  SLAVE  connection  points 
for  implementation  of  the  Phase-Locked  Loop  Board  shown  in  Figure  (6-7).  This  board 
synchronizes  the  switching  frequencies  of  both  Source  Buck  Converter  controllers  so  that 
they  can  operate  in  parallel.  Without  proper  synchronization,  load  sharing  between 
paralleled  units  cannot  be  achieved  due  to  inductor  current  oscillations.  As  shown  in 
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reference  [7],  digital  control  precludes  the  need  for  this  circuit  since  frequency  can  be 
directly  set  by  timer  circuits. 

Implementation  of  the  phase-locked  loop  requires  one  controller  to  be 
designated  the  MASTER  so  that  it  can  set  the  switching  frequency.  The  other  board,  by 
default,  is  the  SLAVE.  In  actual  operation,  the  phase-locked  loop  generates  a  DC  current 
depending  on  the  difference  in  frequencies  of  the  two  controllers.  This  current  sums  with 
the  set  current  of  the  SLAVE  PWM  chip,  directly  influencing  the  switching  frequency  of 
the  SLAVE  controller.  This  cycle  continues  until  frequency  capture  and  lock  results.  In 
addition  to  the  functions  and  interfaces  mentioned  above,  the  PWM  Stage  also 
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provides  pulse-by-pulse  and  overcurrent  protection  circuitry.  These  circuits  are 
documented  in  the  following  Protection  and  Startup  Circuitry  Stage  discussion. 

e.  Protection  and  Startup  Circuitry  Stage 

The  analog  controller  provides  the  Source  Buck  Converters  with  two 
forms  of  protection.  The  first  is  the  pulse-by-pulse  current  limiting  feature  of  the 
UC3637  chip.  This  circuit  protects  the  IGBT  from  exceeding  its  90A  current  rating. 
Referring  to  Figure  (6-6),  pins  12  and  13  make  up  this  circuit.  If  a  fault  occurs  in  the 
buck  converter,  iL/10  will  rise.  If  its  peak  value  gets  to  -6.75  V,  the  driver  output  of  U3 

goes  low,  and  the  IGBT  remains  open  until  peak  inductor  current  goes  below  the  6.75V 
threshold. 

The  second  protection  circuit  is  overcurrent  time-out  shown  in 
Figure  (6-8).  This  circuit  protects  components  from  thermal  damage  when  io/10  exceeds 
150%  of  its  rated  value  for  more  than  1  second.  Utilizing  the  op  amp  in  U3  (pins  15-17), 
an  integrator  was  designed  such  that  io/10  =  4.5  V  causes  pin  17  to  reach  -10V  in  1 
second.  This  voltage  trips  comparator  U7  high  and  causes  U52  and  U54  to  go  high  (~ 

+V).  As  a  result,  comparator  U8  goes  high  initiating  a  shutdown  signal  to  pin  14  of  U3. 
This  protective  action  has  the  end  effect  of  halting  Source  Buck  Converter  operation. 

The  startup  circuitry  allows  duty  cycle  to  ramp  up  to  its  steady-state  value 
from  zero  initial  conditions.  When  the  controller  is  initially  energized,  U5,  goes  high 
due  to  a  +5V  pulse  generated  by  the  RC  circuit  at  input  pin  1  of  the  OR  gate.  As  a  result, 
U54  and  Ug  go  high  and  Q1  turns  “on”.  This  action  causes  pin  3  of  U6  to  go  low  (~  -V) 
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and  prevents  a  reference  signal  from  being  generated.  The  end  result  is  a  voltage  lockout. 
This  action  prevents  a  duty  cycle  waveform  from  being  generated  instantaneously.  If 
such  an  occurrence  happened  from  zero  initial  conditions,  large  current  and  voltage 
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Figure  6-8,  Analog  Control  Board:  Protection  and  Startup  Circuitry  Stage 


71 


oscillations  would  result  with  possible  damage  to  the  buck  converter.  As  it  stands, 
however,  d(t)=0.  By  resetting  momentary  switch  S2,  Q1  is  turned  “off 5  and  voltage  is 
applied  to  pin  3  of  U6.  Reference  voltage  ramps  up  to  its  final  setpoint  through  the  action 
of  the  RC  time  constant  established  by  R603  and  C603  (t  =  1  sec).  Raising  the  reference 
voltage  in  this  way  causes  duty  cycle  to  ramp  to  its  operating  value  thus  preventing  the 
aforementioned  current  and  voltage  oscillations  on  startup.  After  startup,  the  reference 
voltage  continues  to  be  supplied  through  U6. 


f.  Power  Supply  Stage 

Shown  in  Figure  (6-9),  the  power  supply  takes  1 17  VAC,  60  Hz  receptacle 
power,  and  rectifies  it  into  1 1.5  VDC  for  input  to  the  DATEL  TWR-5/1000-15/200-D12 
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voltage  regulator.  The  output  of  the  DATEL  module  is  +5  VDC  and  +13.5  Vdc  supply 
power  for  the  control  board.  The  +15  VDC  is  not  obtained  due  to  minimal  circuit  loading. 
The  lower  +13.5  Vdc  values  meet  control  board  component  power  supply  requirements. 


4.  IGBT  Driver  Board 


The  IGBT  driver  board  completes  the  feedback  loop.  Input  to  the  board  can  be 


made  using  twisted  pair  (analog  control)  or  optical  interface  (digital  control).  Both 


configurations  are  shown  in  Figure  (6-10).  The  input  driver  signal  is  amplified  to  +  10  V 


by  the  optically  isolated  push-pull  transistor  pair  in  U1 .  Following  amplification,  the 
driver  signal  is  applied  to  the  IGBT  gate.  The  back-to-back  1N4744  zener  diodes  (Vz  =  15 


V)  were  installed  across  the  gate  to  prevent  overvoltage  damage. 
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Figure  6-10,  IGBT  Driver  Board 
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5.  Digital  Control  Board 

Reference  [7]  contains  a  detailed  description  on  the  layout  and  operation  of  the 
Load  Buck  Converter  digital  control  board. 
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VII.  CONCLUSIONS 


This  research  documented  the  design  of  SSCM  units  and  their  controllers  for  the 
PEBB  Network  Simulation  Testbed  being  built  at  the  Naval  Postgraduate  School.  In 
addition,  baseline  simulations  were  developed  to  identify  system  instabilities  and 
interesting  controller  interactions.  Finally,  hardware  validation  studies  were  conducted  to 
confirm  the  PEBB  Network  simulations. 

The  SSCM  design  process  began  with  the  selection  of  components  for  the  buck 
converter  power  sections.  The  initial  selections  were  based  upon  system  specifications, 
theoretical  calculations,  and  physical  limitations.  The  process  continued  with  the 
development  of  state-difference  analog  and  digital  control  algorithms.  Next,  gain 
selection  was  performed  for  the  Source  Buck  Converter  analog  controller  and  the  Load 
Buck  Converter  digital  controller.  Chapter  IV  showed  that  to  achieve  the  desired 
response,  proper  output  capacitor  sizing  and  current  gain  selection  had  to  be  achieved  to 
minimize  unwanted  harmonic  effects  in  the  system  output. 

With  the  proper  gain  values  determined,  detailed  ACSL  simulations  were 
performed  as  described  in  Chapter  V.  The  results  showed  that  the  PEBB  Network  is 
stable  under  the  seven  basic  operational  modes.  In  addition,  hardware  validation  studies 
were  conducted  for  Modes  1, 4  and  5.  The  results  of  this  testing  verified  that  the  ACSL 
simulations  correlated  well  with  hardware  implementation.  From  this  process,  it  was 
shown  that  the  SSCM  units  operated  to  their  design  specifications. 
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One  issue  not  considered  in  the  initial  design  process,  however,  dealt  with 
shielding.  As  shown  in  Chapter  V,  switching  noise  had  a  considerable  effect  on  the 
system.  This  problem  will  only  grow  worse  with  the  addition  of  more  SSCM  and  SSIM 
units.  As  a  result,  the  next  generation  of  SSCMs  need  to  employ  methods  to  eliminate 
switching  noise.  Several  recommendations  include: 

•  soft  switching 

•  placement  of  an  enclosure  about  the  buck  converter  power  sections 

•  relocation  of  the  buck  converter  controllers  to  a  separate  enclosure 

•  use  of  fiberoptic  connections  between  controllers  and  feedback  paths. 

An  additional  system  improvement  deals  with  the  replacement  of  the  Source  Buck  analog 
controllers  with  digital  controllers.  This  step  will  greatly  enhance  the  flexibility  of  the 
PEBB  Network  since  frequency  and  gain  modifications  can  be  performed  in  software  for 
the  digital  controller.  As  a  result,  a  wide  variety  of  system  interactions  can  be 
investigated  on  the  fly  rather  than  having  to  wait  for  resistor  replacements  on  the 

analog  control  boards.  In  addition,  the  digital  control  provides  another  layer  of  noise 
immunity. 

Recommendations  for  future  research  focus  on  further  development  of  the  ACSL 
simulation  models.  Areas  for  investigation  include: 

•  effect  of  buck  converters  operating  at  different  frequencies 

•  implementation  of  an  ARCP  synchronous  reference  frame  current  controller 

•  incorporation  of  Gaussian  white  noise  into  the  sensed  variables  to  model  the 
effects  of  switching  noise. 
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Ideally,  these  detailed  simulations  can  then  be  used  to  direct  hardware  testing  towards 
interesting  PEBB  Network  interactions  and  stability  bounds. 
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A. 


APPENDIX  A.  DATASHEETS 

INTERNATIONAL  RECTIFIER  IRGT1090U06 IGBT 


kitemational 
glS  Rectifier 


PD-9.954B 

IRGTI090U06 


■HALF-BRIDGE’  IGBT  INT-A-PAK 


•Rugged  Design 
•Simple  gate-drive 

•Ultra-fast  operation  up  to  25KHz  hard 
switching,  or  1 OOKHz  resonant 
•Switching-Loss  Rating  includes  all  “tail 
losses 


Description 

IR’s  advanced  IGBT  technology  is  the  key  to 
this  line  of  INT-A-pak  Power  Modules.  The 
efficient  geometry  and  unique  processing  of 
the  IGBT  allow  higher  current  densities  than 
comparable  bipolar  power  module  transistors, 
while  at  the  same  time  requiring  the  simpler 
gate-drive  of  the  familiar  power  MOSFET .  This 
superior  technology  has  now  been  coupled  to 
state  of  the  art  assembly  techniques  to  pro¬ 
duce  a  higher  current  module  that  Is  highly 
suited  to  power  applications  such  as  motor 
drives,  uninterruptible  power  supplies,  weld¬ 
ing,  induction  heating  and  ultrasonics. 


Ultra-fast™  Speed  IGBT 


Half-Bridge 

•  *0^ 

V__  =  600V 

C  4 

CE 

lc=90A 

C  a  . 

Vce(ON)<  3.0V 

"  —2 

Parameter  |  Description 

Value 

Units 

v , ,  ,  I  Continuous  collector  to  emitter  voltage 

600 

V 

|c®  Tc  =  25‘C 

Continuous  collector  current 

90 

A 

jc©Tc=  85*C 

Continuous  collector  current 

60 

lc0Tc  =  100°cl 

Continuous  collector  current 

50 

*LM 

Peak  switching  current 

180 

IFM 

Peak  diode  forward  current  (1)  225 

VGF 

Gate  to  emitter  vottage  - 

±20 

V 

Vicru 

RMS  isolation  voltage,  any  terminal  to  case,  t=  1  min 

2500 

pD  ©  Tc  =  25° C  : 

Power  dissipation 

298 

W 

Tj 

Operating  junction  temperature  range 

-40  to  1 50 

°C 

tstq 

Storage  temperature  range 

-40  to  125 

(1)  Duration  limited  by  max  junction  temperature. 


IRGTI090U06 


Electrical  Characteristics  -Tj  =  25,C,  unless  otherwise  stated 


Description 

Min 

Typ  Max 

Co  fleet  or-to -emitter  breakdown  voltage 

600 

—  _ 

CoJiector-to-emitter  voltage 

— 

EBB 

Diode  forward  voltage  -  maximum 

— 

nn 

BSBi 

Gate  threshold  voltage 

El 

BiES 

Threshold  voltage  temperature  coeff. 

BIB 

Forward  transconductance 

a 

caca 

Collector-to-emitter  leakage  current 

— 

BUI 

Gate-to-emitter  leakage  current 

— 

BB1 

Voe  =  i20V 


Dynamic  Characteristics  -  Tj  =  150°C 


Parameter!  Description 


Tum-on  switching  energy 
T um-off  switching  energy 
Total  switching  energy 


Tum-on  delay  time 
Rise  time 

Turn-off  delay  time 
Fall  time 


Diode  peak  recovery  current 
Diode  recovery  time 
Diode  recovery  charge 


Gate-to-emltter  charge  (turn-on)  150 

Gate-to-collector  charge  (tum-on)  70 

Total  gate  charge  (tum-on)  26 


Input  capacitance 
Output  capacitance 
Reverse  transfer  capacitance 


(1) Includes  tail  fosses 

Thermal  and  Mechanical  Characteristics 


VCC  =  30V 
f  =  1MHz 


Parameter 

Description 

Rdwc  flGBT) 

Thermal  resistance,  junction  to  case,  each  1GBT 

R»uc  (Diode) 

Thermal  resistance ,  junction  to  case ,  each  diode 

Rihes  (Module) 

Thermal  resistance ,  case  to  sink 

Wt 

Weight  of  module 

80 


agggsaS9888R889aa 

Ov)  (jairo  PJ0«P)  *x»uo|J3)M 


mi 


INI 


40  80  120  160  2QD  243 

1  10 

Totel  Goto  Choge  <nO 

VCT.  Catecfor-toErrltler  Mdfcge  (V) 

Fig.  7  -  Typical  Gate  Charge  vs. 

Rg.  8  •  Typical  Capacitance  vs. 

Gate-to-Emltter  Voltage 

Col  tector-to- Emitter  Voltage 

25  SO  75  im  125  ISO 
Tj,  Juririlon  Terrperchj©  C°0 
Rg.  9  -  Typical  Switching  Losses 
vs.  Junction  Temperature 


7^j  -lsyc-j- 


0  20  40  dO  80  ' 

I  c,  CoOector-tc>€rnltter  Curent  (A) 
Pig.  10  *  Typical  Switching  Losses  vs. 
Col  lector-to- Emitter  Current 


i — 


10  30  50  TO  90  110  130  150  170 

R  ^  Goto  RwWcrce  CohrrO 
Ftg.  1 1  -  Typical  Switching  Losses 
vs.  Gate  Resistance 


0  100  23D  300  4X1  500  600  700 

V^.  Macaued  c*  TenrrincIsCV) 

Rg.  12  -  Reverse  Bias  Safe  Operating  Area 
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IRGTI090U06 


Rg.  13  -  Typical  RMS  Output  Current  per  phase  vs.  Frequency  (Synthesized  Sinusoidal  Wave) 


9 
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B. 


UNITRODE  UC3637  PWM  DRIVER  IC 


DESCRIPTION 

The  UC1637  is  a  pulse  width  modulator  circuit  intended  to  be  used  for  a  variety  of  PWM  motor 
drive  and  amplifier  applications  requiring  either  uni-directional  or  bi-directional  drive  circuits. 
When  used  to  replace  conventional  drivers,  this  circuit  can  increase  efficiency  and  reduce 
component  costs  for  many  applications.  All  necessary  circuitry  is  included  to  generate  an  analog 
error  signal  and  modulate  two  bi-directional  pulse  train  outputs  in  proportion  to  the  error  signal 
magnitude  and  polarity. 

This  monolithic  device  contains  a  sawtooth  oscillator,  error  amplifier,  and  two  PWM 
comparators  with  ±100mA  output  stages  as  standard  features.  Protection  circuitry  includes 
under-voltage  lockout,  pulse-by-pulse  current  limiting,  and  a  shutdown  port  with  a  2.5V 
temperature  compensated  threshold. 

The  UC1637  is  characterized  for  operation  over  the  full  military  temperature  range  of  -55°C  to 
+125°C,  while  the  UC2637  and  UC3637  are  characterized  for  -25°C  to  +85°C  and  0°C  to  +70°C, 
respectively. 


ABSOLUTE  MAXIMUM  RATINGS  (Note  1) 


Supply  Voltage  (±Vs) . ±20V 

Output  Current,  Source/Sink  (Pins  4,  7) . . 500mA 

Analog  Inputs  (Pins  1,  2, 3,  8, 9,  10,  1 1  12, 13,  14, 15, 16) . ±Vs 

Error  Amplifier  Output  Current  (Pin  17) . ±20mA 

Oscillator  Charging  Current  (Pin  18) . -2mA 

Power  Dissipation  at  TA  =  25°C  (Note  2) .  1  OOOmW 

Power  Dissipation  at  TC  =  25°C  (Note  2) .  2000mW 

Storage  Temperature  Range . -65°C  to  +150°C 

Lead  Temperature  (Soldering,  10  Seconds) . +300°C 


Note  1:  Currents  are  positive  into,  negative  out  of  the  specified  terminal. 

Note  2:  Consult  Packaging  Section  of  Databook  for  thermal  limitations  and  considerations  of package. 

BLOCK  DIAGRAM 
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CONNECTION  DIAGRAM 


ELECTRICAL  CHARACTERISTICS: 

Unless  otherwise  stated,  these  specifications  apply  for  TA  =  -55°C  to  +125°C  for  the  UC1637-  -25°r  to 
+85°C  ^  the  UC2637;  and  0°C  to  +70°C  for  the  UC3637;  +VS  =+15V,  -VS  =  15V  Sra  =  5V  VTH 
=  -5V,  RT  =  16.7kW,  CT  =  1500pF,  TA=TJ.  ’  ™  5V>-VTH 
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ELECTRICAL  CHARACTERISTICS: 

Unless  otherwise  stated,  these  specifications  apply  for  TA  =  -5 5°C  to  +125°C  for  the  UC1637;  -25°C  to 
+85°C  for  the  UC2637;  and  0°C  to  +70°C  for  the  UC3637:  VS  =+15V,  -VS  =  -  15V,  +VTH  =  5V,  -VTH  = 
-5V,  RT  -  16.7kW,  CT  =  1500pF,  TA=TJ. 


tssTcoNornoNs 

!  UC1637/UC2637 

UC3S37 

UNITS 

Eaizaiag 

fTlTBKjTji;T:V< 

|EnorAmpifter{Cominu«<0  1 

Ouipui  Sink  Current 

VpW  17*0V 

1  5 

-50 

B2^3B 

Output  Source  Current 

Vpw  17*»0V 

mm 

L_ 

5 

11 

_ i 

Hegfc  Level  Output  Vcltasp 

mm 

BUM 

t” 

13- 

13.6 

_ 

mm 

Low  Level  Output  Voltage 

□ 

BB3 

KM 

EEH 

KM 

mm 

I  PWM  Comparators  I 

zn 

VCM»6V 

20 

20 

_ 

KM 

Input  Bias  Current 

Vcm^CV 

mm 

KM 

wm 

■a 

ISM 

Input  bfystaress 

Vcm  m  GV 

mm 

mm 

KM 

Common  Mode  range 

Vs  »±SVtoi20V 

ESB 

BBS! 

ESS 

n 

|  Current  Limit  | 

V0M  =  QV.Tj=25‘C 

1131 

ESI 

EH 

KM 

■ 

E9 

EH 

RMIMl 

mem 

1t>M 

EH 

m A 

I  Common  Mode  Range 

Vfis  ±2-5Vto±20V 

E3 

KM 

RBI 

V 

|  Shutdown  ! 

1  Shutdown  Threshold 

(Ntote4) 

Kiel 

ej 

mm 

*£.5 

mm 

i£j 

B.l  li  1 1— 

40 

KM; 

I  Input  Bias  Current 

VPWH  =  ^VStO‘VS 

KM 

-0.5 

iiA 

|  Under-Voltage  Lockout  | 

I  Start  Threshold 

IIISXHHHHHHi 

K39 

4,15 

El 

mm 

0,25 

mm 

KMl 

Total  Standby  Current  | 

Output  Low  Level 

l$tNK=20mA 

mu 

mm 

dEli 

■9 

V 

(sink  =  ioomA 

mm 

BRI 

KM 

Output  Level 

ISCU^  =  20mA 

mm 

mm 

135 

V 

!  source  s  100mA 

mm 

B«U 

mm 

13.5 

Rise  Time 

(Note  3)  a  » Inf,  Tj«  ZS'C 

100  1 

600 

KF71 

ro 

FalTime 

(Nets  3)0.=  Inf,  Tj=25'C 

KB1 

BSB1 

El 

mm 

Note  3:  These  parameters,  although  guaranteed  over  the  recommended  operating  conditions ,  are  not  100%  tested 
in  production. 

Note  4:  Parameter  measured  with  respect  to  +  VS  (Pin  6). 

Note  5:  Parameter  measured  at  +VS  ( Pin  6)  with  respect  to  -VS  (Pin  5). 

Note  6:  RT  and  CT  referenced  to  Ground. 


FUNCTIONAL  DESCRIPTION 

Following  is  a  description  of  each  of  the  functional  blocks  shown  in  the  Block  Diagram. 
Oscillator 

The  oscillator  consists  of  two  comparators,  a  charging  and  discharging  current  source,  a  current 
source  set  terminal,  1SET  and  a  flip-flop.  The  upper  and  lower  threshold  of  the  oscillator 
waveform  is  set  externally  by  applying  a  voltage  at  pins  +VTH  and  -VTH  respectively.  The 
+VTH  terminal  voltage  is  buffered  internally  and  also  applied  to  the  ISET  terminal  to  develop 
the  capacitor  charging  current  through  RT.  If  RT  is  referenced  to  -VS  as  shown  in  Figure  1,  both 
the  threshold  voltage  and  charging  current  will  vary  proportionally  to  the  supply  differential,  and 
the  oscillator  frequency  will  remain  constant.  The  triangle  waveform  oscillators  frequency  and 
voltage  amplitude  is  determined  by  the  external  components  using  the  formulas  given  in  Fig.  1. 
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FigKir*  1.  Oseiitetor  Setup 


PWM  Comparators 

Two  comparators  are  provided  to  perform  pulse  width  modulation  for  each  of  the  output  drivers. 
Inputs  are  uncommitted  to  allow  maximum  flexibility.  The  pulse  width  of  the  outputs  A  and  B  is 
a  function  of  the  sign  and  amplitude  of  the  error  signal.  A  negative  signal  at  Pin  10  and  8  will 
lengthen  the  high  state  of  output  A  and  shorten  the  high  state  of  output  B.  Likewise,  a  positive 
error  signal  reverses  the  procedure.  Typically,  the  oscillator  waveform  is  compared  against  the 
summation  of  the  error  signal  and  the  level  set  on  Pin  9  and  11. 

MODULATION  SCHEMES 

Case  A  Zero  Deadtime  (Equal  voltage  on  Pin  9  and  Pin  1 1) 

In  this  configuration,  maximum  holding  torque  or  stiffness  and  position  accuracy  is  achieved 
However,  the  power  input  into  the  motor  is  increased.  Figure  3A  shows  this  configuration. 

Case  B  Small  Deadtime  (Voltage  on  Pin  9  >  Pin  1 1) 

A  small  differential  voltage  between  Pin  9  and  1 1  provides  the  necessary  time  delay  to  reduce 
the  chances  of  momentary  short  circuit  in  the  output  stage  during  transitions,  especially  where 
power-amplifiers  are  used.  Refer  to  Figure  3B. 

Case  C  Increased  Deadtime  and  Deadband  Mode  (Voltage  on  Pin  9  >  Pin  1 1) 

With  the  reduction  of  stiffness  and  position  accuracy,  the  power  input  into  the  motor  around  the 
null  point  of  the  servo  loop  can  be  reduced  or  eliminated  by  widening  the  window  of  the 
comparator  circuit  to  a  degree  of  acceptance.  Where  position  accuracy  and  mechanical  stiffness 
is  unimportant,  deadband  operation  can  be  used.  This  is  shown  in  Figure  3C. 
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Output  Drivers 

Each  output  driver  is  capable  of  both  sourcing  and  sinking  100mA  steady  state  and  up  to  500mA 
on  a  pulsed  basis  for  rapid  switching  of  either  POWERFET  or  bipolar  transistors.  Output  levels 
are  typically  -VS  +  0.2V  @50mA  low  level  and  +VS  -  2.0V  @50mA  high  level. 

Error  Amplifier 

The  error  amplifier  consists  of  a  high  slew  rate  (15V/ms)  op-amp  with  a  typical  1MHz 
bandwidth  and  low  output  impedance.  Depending  on  the  ±VS  supply  voltage,  the  common  mode 
input  range  and  the  voltage  output  swing  is  within  2V  of  the  VS  supply. 

Under-Voltage  Lockout 

An  under-voltage  lockout  circuit  holds  the  outputs  in  the  low  state  until  a  minimum  of  4V  is 
reached.  At  this  point,  all  internal  circuitry  is  functional  and  the  output  drivers  are  enabled.  If 
external  circuitiy  requires  a  higher  starting  voltage,  an  over-riding  voltage  can  be  programmed 
through  the  shutdown  terminal  as  shown  in  Figure  4. 


Shutdown  Comparator 

The  shutdown  terminal  may  be  used  for  implementing  various  shutdown  and  protection  schemes. 
By  pulling  the  terminal  more  than  2.5V  below  VIN,  the  output  drivers  will  be  enabled.  This  can 
be  realized  using  an  open  collector  gate  or  NPN  transistor  biased  to  either  ground  or  the  negative 
supply.  Since  the  threshold  is  temperature  stabilized,  the  comparator  can  be  used  as  an  accurate 
low  voltage  lockout  (Figure  4)  and/or  delayed  start  as  in  Figure  5.  In  the  shutdown  mode  the 
outputs  are  held  in  the  low  state. 
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Current  Limit 

A  latched  current  limit  amplifier  with  an  internal  200mV  offset  is  provided  to  allow  pulse-by- 
pulse  current  limiting.  Differential  inputs  will  accept  common  mode  signals  from  -VS  to  within 

3V  ofthe  +vs  suPply  while  providing  excellent  noise  rejection.  Figure  6  shows  a  typical  current 
sense  circuit. 


FlguwS.  Cterert  Limit  Sensing 


90 


APPENDIX  B.  SIMULINK  BLOCK  DIAGRAMS  AND  MATLAB  CODE 


A.  SOURCE  BUCK  CONVERTER 

The  Matlab  m-files  and  Simulink  block  diagram  used  to  generate  first-cut  gain 
selections  for  the  Source  Buck  Converter  are  provided  below. 


1.  Matlab  M-file  buckcons.m 


%************************ ***************** ******* 
%  buckcons.m 
% 

%  Specifications: 

% 

%  90  A  IGBPs 

%  Analog  controller  patterned  after  SSCM. 

% 

%  Want  continuous  operations  between  10%  and 
%  100%  loading. 

% 

%  Simulation  sequence  cycles  through 
%  startup  ->  10%  loading  ->  full  loading. 

% 

%  Written  by:  Mike  Badorf 
%  Last  mod:  2JUN97 

o ^* *********************************************** 


%******************** **************** ************ 
%  Constant  Definitions 

%*****************************  **********  ********* 


E 

=  400.0; 

%  input  voltage  of  buck 

Vcd 

=  305.0; 

%  desired  output  voltage  of  buck:  10% 

Pout 

=  9e3; 

%  output  power 

Cf 

=  400e-6; 

%  chosen  based  on  fres  «  fline  &  parts 

Lfilt 

=  0.425e-3; 

%  actual  filter  inductance  value 

Cfilt 

=  2000e-6; 

%  actual  filter  capacitance  value 

rxfull 

=  9.68; 

%  low  R  value  -  full  load  current 

rx_ten 

=  96.8; 

%  high  R  value  - 10%  load  current 

rxstart 

=  rx_ten; 

%  starting  resistance 

d 

=  Vcd/E; 

%  duty  cycle 

time 

=  360.0; 

%  360  Hz  ripple  from  6  pulse  rectifier 

f 

=  20e3; 

%  switching  frequency 

T 

=  l/f; 

%  switching  period 

pw 

=  d*T; 

%  pulse  width  for  switching  noise  SS  model 

L 

=  0.760e-3; 

%  actual  inductance  value  for  buck  1 

<^******* *************  ********************  ******** 
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%  Determination  of  Lcrit:  low  load  R  determines 

%  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  sf:  *  *  *  *  %  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  * 

R  =  rx  ten; 

Lcrit  =  (T*R/2)*(l-d);  %  critical  inductance 


°/o  *  *  *  *  *  sfc  *  *  *  *  *  *  afe  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  4c  *  *  sjc  * 

%  Determination  of  deltaj,  Imax,  Imin,  ILavg 
%  and  I  Joad  at  1 0%  load  (11). 

%*** * * ************* **** ***** ** * * # * * * * * * * ** * * * * * * * 


deltaj  =  ((E  -  Vcd)/L)*d*T; 

Imax  11  =  d*E*(l/R  +  (l-d)*T/(2*L)); 
Imin  11  =  d*E*(l/R  -  (l-d)*T/(2*L)); 
ILavg  11  =  (Imax Jl+IminJl)/2; 
IJoadJl  =  Vcd/R; 


%  II  max  L  current 
%  11  min  L  current 
%  11  avg  L  current 
%  11  current 


%  Determination  of  Imax,  Imin,  ILavg 
%  and  I  load  at  lull  load  (fl). 


R  =  rx  full; 

Imax  =  d*E*(l/R  +  (l-d)*T/(2*L)); 
Imin  =  d*E*(l/R  -  (l-d)*T/(2*L)); 
ILavg  =  (Imax+Imin)/2; 

I  Joad  =  Vcd/R; 


%  fl  max  L  current 
%  fl  min  L  current 

%  fl  avg  L  current 
%  fl  current 


%  Determination  of  resonant  peak  for  L  /  CfLPF. 
fres  =  l/(2*pi*sqrt(L*Cf)); 


%  Determination  of  resonant  peak  for 
%  input  Lfilt  /  Cfilt  LPF. 

%♦**** ***************************^t ******  ******** 


fres_filt  =  l/(2*pi*sqrt(Lfilt*Cfilt)); 


I 
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2. 


Matlab  M-file  buckstrt.m 


0^  afe  *  afc  sfc  $  sjc  *  *  *  :fc  jjc  sfc  $  *  *  gfe  $  j|c  $  $  jjc  $  afc  $  $  *  $  jjt  $  s|c  jf:  +  *  ajc  a|e  sjc  *  *  *  jfc  *  $  *  afe  *  jjc  J|e  * 


%  buckstrt.m 


% 

%  Provides  startup  parameters  for  9  kW,  20  kHz 
%  buck  chopper. 

% 


%  Calls  buckcons.m,  a  constant  file  developed 
%  for  above  buck  chopper. 

% 

%  Written  by:  Mike  Badorf 
%  Last  mod:  2JUN97 


%**** ********************************** ********** 


clear  %  clear  MATLAB  memory 
buckcons  %  establish  buck  constants 


o/0***** ********************** ********************* 

%  Bessel  Pole  specifications 

( ^** ************************* ****** *************** 
wo  =  3250; 

poles  =  wo*[  -0.942  -0.74550+0.71 12*i  -0.7455-0.71 12*i]; 

o^*^**ic****** **************** ********************* 

%  desired  characteristic  polynomial  used  to 
%  determine  gains 

%**  *  **  *  ***  *  *  ***  **  **  ****  *  ****  *  *  *  *  **  *  *  *  *  ****  **  *  **** 
S  =  poly(poles); 

o^** ******* ***** ******* ********* ******* *** *** ***** 

%  Gain  calculations:  Calculated  using  low 
%  starting  resistance.  Gains  are  utilized  as 
%  first  run  starting  point.  Entry  of  manual 
%  gains  can  follow  to  give  best  overall 
%  response  over  all  operating  conditions.  Use 
%  board  resistors  in  practice  to  set  gains. 
%************************************************ 
R  =  rxstart; 

hi  =  ((L*R*Cf*S(2))-L)/(E*R*Cf) 

hn  =  (L*Cf*S(4))/E  +  hi*L/2 

hv  =  ((L*R*Cf*S(3»-(E*hi)-R)/(E*R) 

0,/ ************************************************ 

%  Function  to  verify  poles  obtained  with  gains 
%  match  desired  poles. 

%* ************************************ *********** 
Am  =  [  0  -1/L  E/L 

1/Cf  -l/(R*Cf)  0 

-hv/Cf  (h  i/L)+(h  v/(R*  Cf))-hn  -hi*E/L  ]; 


95 


closedsyspole  -  eig(Am)  %  check  for  correct  poles 

0/^444444444444444444444444444444444444444444444444 
%  SIMULINK  no  load  parameters  (are  in  addition 
%  to  constants  determined  in  buckcons.m) 

%444444444444444444444444444444444444444444444444 


%  ABCD  state  space  parameters 
A  =  [  -l/(R*Cf)  1/Cf 
-1/L  0  ]; 

B  =  [  0 
E/L]; 

C  =  [  1  0 

0  1]; 

D  =  [  0 

0]; 


<^444444444444444444444444444444444444444444444444 

%  Update  initial  conditions  for  integrators 
0^444444444444444444444444444444444444444444444444 


=  304; 

%  initial  Vc  integrator  value 

=  Vciv/R; 

%  initial  hn  integrator  value 

/  =  0.0; 

%  initial  hn  integrator  value 

=  304/E; 

%  initial  duty  cycle  value 

=  0.0; 

%  simulation  start  time 

=  0.012; 

%  simulation  stop  time 
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3. 


Matlab  M-file  buckfull.m 


<^***************jm  ************************* ****** 


%  buckfulkm 


% 

%  Provides  full  load  parameters  for  9  kW, 
%  20  kHz  buck  chopper. 

% 


%  Written  by:  Mike  Badorf 


%  Last  mod:  2JUN97 

0^**** ************************************ ******** 


R  =  rx  full;  %  full  load  resistance  value 


****************************** **************** 


%  Function  to  verify  poles  obtained  with  gains 
%  match  desired  poles. 

%*** ********************************** *********** 


Am-[0  -1/L  E/L 

1/Cf  -l/(R*Cf)  0 

-hv/Cf  (hi/L)+(hv/(R*Cf))-hn  -hi*  E/L  ]; 


syspole  =  eig(Am);  %  check  for  correct  poles 

%************************************************ 
%  SIMULINK  no  load  parameters  (are  in  addition 
%  to  constants  determined  in  buckcons.m) 

<^********* *************************************** 


%  ABCD  state  space  parameters 
A  =  [ -l/(R*Cf)  1/Cf ;  -1/L  0  ];B-[0;  E/L]; 

C  =  [  1  0;  0  1 1; D  =  [  0;  0]; 


%  Data  passdown  before  buckfull.m  simulation 
%  begins. 

o/0** t  ***************************************  ******* 


Vc_response  =  Vc; 
iL_response  -  iL; 
duty_response  =  duty; 
timer  =  time; 


%*********************************** ************* 


%  Update  initial  conditions  for  integrators  in 
%  buckfull.m  simulation. 

%  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  **  *  *  *  *  *  *  *  *  *  *  **  *  *  *  *  *  *  *  * 

Vciv  =  Vc(length(time));  %  Vc  initial  value  for  buckfull.m 
iLiv  =  iL(length(time));  %  iL  initial  value  for  buckfull.m 
dutyiv  =  duty(length(time));  %  duty  initial  value  for  buckfull.m 
tstart  =  tstop;  %  sim.  start  time  for  buckfull.m 

tstop  =  tstop  +  0.012;  %  simulation  stop  time  for  buckfull.m 
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4. 


Matlab  M-file  buckten.m 


%*  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  afc  ****  **  **  *  *  **  *  *  *  *  *  *  *  *  *  *  ***  * 

%  buckten.m 
% 

%  Provides  10%  load  parameters  for  9  kW,  20  kHz 
%  buck  chopper. 

% 

%  Written  by:  Mike  Badorf 
%  Last  mod:  2JUN97 

R  =  rx__ten;  %  10%  load  resistance  value 

%  ***********  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  4c  *  9|c  *  *  *  *  *  *  *  * 

%  Function  to  verify  poles  obtained  with  gains 
%  match  desired  poles. 

%  ******  *  ****  *  ***  *  *  *  *  *  *  *  *  *  *  %  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  * 
Am  =  [  0  -1/L  E/L 

1/Cf  -l/(R*Cf)  0 

-hv/Cf  (hi/L)+(hv/ (R*  Cf))-hn  -hi*E/L]; 

syspole  =  eig(Am);  %  check  for  correct  poles 

%  *  ******  *  *******  *  *  *  *  4c  *  *  *  *  *  *  *  *  3fe  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  * 

%  SIMULINK  no  load  parameters  (are  in  addition 
%  to  constants  determined  in  buckcons.m) 

%**  *  **  *  *  ***  *  ***  *  *  *  **  *  *  ***  *  *  *  *  *  **  *  *  *  *  *  *  ***  *  *  *  *  **  *  * 


%  ABCD  state  space  parameters 
A  =  [  -l/(R*Cf)  1/Cf;  -1/L  0  ];B=[0;  E/L]: 

C  =  [10;  0  1  ];  D  =  [  0;  0]; 


%  *******  *  ***  *  ****  *  *  *  *  *  *  *  *  *  %  4c  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  * 

%  Data  passdown  before  buck__ten.m  simulation 
%  begins. 

%  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  %  *  *  * 

Vcjesponse  =  [  Vcjesponse;  Vc]; 

iL_response  =  [  iL  response;  iL]; 

duty_response  =  [  dutyresponse;  duty]; 

timer  =  [  timer;  time]; 


o/o***********************************^^^^* 

%  Update  initial  conditions  for  integrators  for 
%  buckten.m  simulation. 

V civ  =  Vc(length(time));  %  Vc  initial  value  for  buckten.m 
iLiv  =  iL(length(time));  %  iL  initial  value  for  buckten.m 
dutyiv  =  duty(length(time»;  %  duty  initial  value  for  buckten.m 
tstart  =  tstop;  %  simulation  start  time  for  buckten.m 

tstop  =  tstop  +  0.0 1 2;  %  simulation  stop  time  for  buckten.m 
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5. 


Matlab  M-file  buckplot.m 


0/0  ************************************************ 


%  buckplot.m 
% 

%  Plot  startup,  change  to  100%  loading,  & 

%  change  to  10%  loading  for  9  kW,  20  kHz  buck 
%  chopper. 

% 


%  Written  by:  MikeBadorf 
%  Last  mod:  2 1 SEP96 


<^/*  **********************  ************ ************* 


Vc__response  =  [  Vcjresponse;  Vc  ]; 
iL_response  =  [  iL_response;  iL  ]; 
duty_response  =  [  duty_response;  duty]; 
timer  =  1 000*  [  timer;  time  ] ; 

datastep  =1:1  :length(timer);  %  Allows  plotting  fewer  points 

%  than  generated. 


figure(l) 

subplot(3,l,l),  plot(timer(datastep),  Vc_response(datastep),'b'),  grid 
title('Vc  Response’), ylabel('V c  (volts)') 

subplot(3,l,2),  plot(timer(datastep),  iL_response(datastep),'b'),  grid 
title('iL  Response'), ylabel('iL  (amps)') 

subplot(3,l,3),  plot(timer(datastep),  10*duty_response(datastep),'b'),  grid 
title('d(t)  Response’), xlabel('time  (msec)’),ylabel('10*d(t)’) 
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6.  Simulink  Block  Diagram 


Figure  B-l,  Source  Buck  Converter  Simulink  Block  Diagram 
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B.  LOAD  BUCK  CONVERTER 

The  Matlab  m-file,  buckcons.m,  and  the  Simulink  block  diagram  for  the  Load 
Buck  Converter  are  provided  below.  The  remaining  m-files  are  provided  above  since  the 
transient  response  generation  is  the  same  for  both  the  Source  and  Load  Buck  Converters. 


1.  Matlab  M-file  buckcons.m 


%*************** ********************* ************ 


%  buckcons.m 


% 


%  Determining  Buck  Circuit  Parameters 
% 

%  Specifications: 

% 


%  90  A  IGBTs 

%  Analog  controller  patterned  after  SSCM. 

% 

%  Want  continuous  operations  between  10%  and 
%  100%  loading. 

% 

%  Simulation  sequence  cycles  through 
%  startup  ->  10%  loading  ->  full  loading. 

% 


%  Written  by:  Mike  Badorf 
%  Last  mod:  21SEP96 

O^***** **** * ******* ******************************* 


%************************************************ 
%  Constant  Definitions 

0/^*  *********************************************** 


E 

=  300.0; 

%  input  voltage  of  buck 

Vcd 

=  208.0; 

%  desired  output  voltage  of  buck 

Pout 

=  3e3; 

%  output  power 

Cf 

=  400e-6; 

%  chosen  based  on  fres  «  fline  &  parts 

Lfilt 

=  0.425e-3; 

%  actual  filter  inductance  value 

Cfilt 

=  460e-6; 

%  actual  filter  capacitance  value 

rx_full 

=  VcdA2/Pout; 

%  low  R  value  -  full  load  current 

rxten 

=  10*rx_full; 

%  high  R  value  -  10%  load  current 

rxstart 

=  rx  ten; 

%  starting  resistance 

d 

=  Vcd/E; 

%  duty  cycle 

fline 

=  360.0; 

%  360  Hz  ripple  from  6  pulse  rectifier 

f 

=  20e3; 

%  switching  frequency 

T 

=  l/f; 

%  switching  period 

pw 

=  d*T; 

%  pulse  width  for  switching  noise  SS  model 

L 

=  1.3e-3; 

%  actual  inductance  value  for  buck 
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%  Determination  of  Lcrit:  low  load  R  determines 
%*****************»****** ****** ****************** 
R  =  rxten; 

Lcrit  =  (T*R/2)*(l-d);  %  critical  inductance 


%  Determination  of  delta_I,  Imax,  Imin,  ILavg 
%  and  I_load  at  10%  load  (11). 


deltaj  =  ((E  -  Vcd)/L)*d*T; 

Imaxll  =  d*E*(l/R  +  (l-d)*T/(2*L)); 
Iminll  =  d*E*(l/R  -  (l-d)*T/(2*L)); 
ILavgll  =  (Imax_ll+Imin_ll)/2; 

I JoadJl  =  Vcd/R; 


%  11  max  L  current 
%  11  min  L  current 
%  11  avg  L  current 
%  11  current 


%  Determination  of  Imax,  Imin,  ILavg 
%  and  I_load  at  full  load  (fl). 


R  =  rx_full; 

Imax  =  d*E*(l/R  +  (l-d)*T/(2*L)); 
Imin  =  d*E*(l/R  -  (l-d)*T/(2*L)); 
ILavg  =  (Imax+Imin)/2; 

IJoad  =  Vcd/R; 


%  fl  max  L  current 
%  fl  min  L  current 
%  fl  avg  L  current 
%  fl  current 


%  Determination  of  resonant  peak  for  L  /  CfLPF. 
fres  =  l/(2*pi*sqrt(L*Cf)); 


%*************************!„+**+#Nt<1#i)[+#+<t<[Nt+](.i(ci).j|ti)[!|i 
%  Determination  of  resonant  peak  for 
%  input  Lfilt  /  Cfilt  LPF. 

%******♦*******  ********111******#^*#***^^^^^^;,,^^^^^ 


fres_filt  =  l/(2*pi*sqrt(Lfilt*Cfilt)); 
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2. 


Simulink  Block  Diagram 


APPENDIX  C.  ACSL  CODE 


A.  SOURCE  BUCK  CONVERTER  GAIN  SELECTION  CODE 

The  computer  code  for  this  section  is  comprised  of  the  main  ACSL  program,  two 
supporting  macro  files  and  the  command  file  used  to  generate  transient  load  studies. 


1.  Main  Program 


!  sourcebuckxsl 

I 

!  This  program  simulates  the  operation  of  a  9  kW  source  buck  converter  unit  operating  at 
!  20  kHz  assuming  ideal  components  and  no  line  losses.  The  converter  has  been  designed 
!  to  operate  as  follows: 

!  400  Vdc  input  voltage 

!  300  Vdc  nominal  output  voltage 

!  continuous  ind.  current  from  10%  ~>  100%  load 

I 

!  The  controller  for  this  simulation  makes  use  of  the  analog  control  algorithm  derived 
!  in  Chapter  III.  The  file,  sourcebuckxmd,  contains  the  runtime  commands  utilized  to 
!  exercise  this  model. 

i 

!  Written  by:  Michael  Badorf 
!  Last  modified:  02JUN97 

INCLUDE  'buck.mac' 

INCLUDE  ’analog  controlmac’ 

PROGRAM  sourcebuck 


INITIAL 

MAXTERVAL  maxt  =  1.0e-6 
MINTERVAL  mint  =  1.0e-8 


CINTERVAL  cint  =  1.0e-6 
ALGORITHM  ialg  =  5 

NSTEPS  nstp  =  1 
CONSTANT  tstop  =  1.0 
CONSTANT  Vramppk  =  10.0 
CONSTANT  swfreql  =20000.0 
Trampl  =  1.0/swfreql 
END  !  "of  initial" 


!  "maximum  integration  step  size" 

!  "minimum  integration  step  size" 

!"for  variable  step  algorithms  " 

!"data  communication  interval" 

{"integration  algorithm" 

!  "4— R.K.  2nd,  5~R.K.4th" 

!  "stop  point  for  integration  " 

!"pk  value  of  ramp  modulating  waveform  " 
{"switching  frequency  for  buck  1  " 

{"switching  period  for  buck  1  " 
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DYNAMIC 


TERMT  (t  .GE.  (tstop-0.5*cint)> 

DERIVATIVE 

Vincapl  =  400.0  !"define  a  hard  input  voltage  of 400  Vdc" 


!  source  buck  converter  unit  1 

|*****.M******#*>M*********iM********:M******iM**#******:M******************** 
tpwml  =  MOD(t, Tramp  1)  ["time  elapsed  in  switching  period  for  buck  1  " 

Vrampl  =  Vramppk*tpwml/Trampl  ! "establish  ramp  modulating  waveform  for  buck  1" 

!  "implement  source  buck  converter  unit  1  macro" 

buck(  1  ,V incap  1  ,iout  1  ,isw  1  ,iL  1  ,Vc, Vramp  1  /'Lbuck  1 =760.0e-6", " Vdidrop  1  =2  0"  & 
"Vswdropl=2.0") 


.'"implement  source  buck  converter  unit  1  multiloop  control  macro 
analog_control(  1 , Vincap  1 , V c,iout  1  ,iL  1 ,  "RCstart  1 =0.6", "hi  1  =0.0 1 5  ","hv  1  =0  0 1 7"  & 
"hnl=26.09") 


CONSTANT  Cbuckl  =  400.0e-6 
CONSTANT  Vcic  =  0.0 
pVc  =  (iLl  -  ioutl)/Cbuckl 
Vc  =  INTEG(pVc,  Vcic) 


("output  capcitance  " 
("output  voltage  IC  " 
("output  voltage  state" 
("output  voltage  " 


CONSTANT  Rout  -  96.8  ("source  buck  converter  startup  load" 

("10%  load  value  " 

ioutl  =  Vc/Rout  ("load  current" 

!  end  source  buck  converter  unit  1 


END  ("of  derivative" 
END  ("of dynamic" 
END  ("of program" 
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2. 


Buck  Converter  Macro 


f  +  ^jjc  +  +  ***3fe**jfC**j|t*5|C******3|C**)|C*******3jc3|c********J|C***>fC******5|C**a|c3»t***5jC***3|t********3|t*3jt5fC***** 

!  buck.mac 

I 

!  This  macro  governs  the  conduction  status  of  a  buck  converter  switch  and  uses  this 
!  condition  to  develop  inductor  current  utilizing  the  inductor  equation  of  state. 

!  Switch  and  diode  losses  are  modeled  as  voltage  drops.  All  other  components  are  ideal. 

!  The  output  inductor  current  is  used  in  the  main  program  to  determine  output  capacitor 
!  voltage.  Although  this  step  could  be  done  in  the  macro,  the  current  arrangement 
!  allows  modeling  more  complex  PEBB  network  configurations. 

!  Written  by:  Michael  Badorf 
!  Last  modified:  02JUN97 

i************************************************************************************* 

MACRO  buck(z,vin,iout,isw,iL,Vc, rampwav, pLbuck,pVdidrop,pVswdrop) 

INITIAL 

CONSTANT  pLbuck 
CONSTANT  pVdidrop 
CONSTANT  pVswdrop 

!  "assign  the  initial  inductor  current" 

CONSTANT  iLic&z  =  0.0 

!  "define  and  initialize  the  switch  conduction  status" 

LOGICAL  SWl&z 

SWl&z  =  .true.  !"T  si  conducts,  F  si  blocks" 


END 


["determine  conduction  status  of  SW1" 
PROCEDURAL(S  W 1  &z,isw=rampwav,duty&z,iL) 

["charge  inductor  when  duty  GT  reframp" 
IF  (duty&z  .GT.  rampwav)  THEN 

SWl&z  =  .true, 
isw  =  iL 


ELSE 

SWl&z  =  .false, 
isw  =  0.0 

ENDIF  ["end  if  (duty&z  .GT.  rampwav) 
END  !"of  procedural" 
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!"Given  conduction  status  of  SW1 ,  and  whether  continuous  or  not,  " 

!"  determine  derivative  of  iL  " 

PROCEDURAL  (piL&z=SWl&z,vin,Vc,Vdidrop&z,iL,Vswdrop&z,Lbuck&z) 

" — if  SW1  gated,  assume  can  conduct  pos.  current" 

IF  (SWl&z)  THEN 

contmode&z  =  .true. 

piL&z  =  (vin-Vswdrop&z-Vc)/Lbuck&z 


ELSE 


IF  (iL  .GT.  0.0)  THEN 

!  "if  IGBT  not  gated  but  iL  greater  than  zero  " 
!"the  diode  conducts  and  the  inductor  discharges" 
piL&z  =  (-Vdidrop&z-Vc)/Lbuck&z 


ELSE 

!  "if  IGBT  not  gated  and  iL  less  than  or  equal" 
!"to  zero,  the  diode  is  off  and  the  inductor  " 
!  "current  and  derivative  are  zero  " 

piL&z  =  0.0 

ENDIF  !"end  if  (iL  .GT.  0.0)" 

ENDIF  !"end  if  (SWl&z)" 

END  !  "of  procedural" 

("integrate  the  inductor  state  variable" 

iL  =  BOUND(0.0,  100.0,  LIMINT(piL&z,  iLic&z,  0.0,  100.0)) 

MACRO  END 
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3.  Analog  Controller  Macro 


(* ********** ****************************************** ************************** ****** 
!  analog_control.mac 

I 

!  This  macro  governs  the  analog  control  algorithm  derived  in  Chapter  III.  In  addition 
!  overcurrent  protection  is  provided  to  place  the  controller  in  pulse-by-pulse 
!  operations  if  inductor  current  exceeds  67.5  amps. 

i 

!  Written  by:  Michael  Badorf 
!  Last  modified:  02JUN97 

MACRO  analog_control(z,Vin,Vout,iout,iL,pRCstart,phi,phv,phn) 

INITIAL 

CONSTANT  pRCstart 
CONSTANT  phi 
CONSTANT  phv 
CONSTANT  phn 

!  "assign  internal  macro  constants" 

CONSTANT  Vthouseref&z  =  305.0 
CONSTANT  houseslope&z  =  3.0 
CONSTANT  Vthouseic&z  =  0.0 
CONSTANT  Vxic&z  =  0.0 
CONSTANT  iLlimit&z  =67.5 
CONSTANT  dutylow&z  =0.0 
CONSTANT  minduty&z  =  0.0 
CONSTANT  maxduty&z  =  10.0 
CONSTANT  scale&z  =  10.0 

END 


!"top  of  housecurve" 
!"housecurve  slope  " 


!"max.  ind.  current  before  prot.  action" 


pVthouse&z  =  (Vthouseref&z- Vthouse&z)/RCstart&z  !  "reference  voltage  state 

Vthouse&z  =  INTEG(pVthouse&z,  Vthouseic&z)  ["reference  voltage 

vcref&z  =  Vthouse&z  -  iout/houseslope&z 

pVx&z  =  Vout  -  vcref&z  ["voltage  state  difference  " 

Vx&z  =  INTEG(pVx&z,  Vxic&z)  ["integral  voltage  state  diff." 

!  "control  algorithm  for  source  buck  converter  controller" 

dutyML&z  =  -hi&z*(iL-iout)-hv&z*pVx&z-hn&z*Vx&z+Vthouse&z/Vin 

["overcurrent  protection  for  buck  converter  controller " 

!"if  buck  converter  1  is  operating,  implement  overcurrent  protection  " 

IF  (iL  .GT.  iLlimit&z)  THEN 

duty&z  =  dutylow&z  !  "overcurrent" 

ELSE 

duty&z  =  BOUND(minduty&z,  maxduty&z,  scale&z*  dutyML&z) 

ENDIF  ["end  if(iL.GT.  iLlimit&z)" 

MACRO  END 
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4. 


Command  File 


******************************t^t*^^Mn,imf^t***************t**************^ii^**^^ii,it:t:,nl 

sourcebuck.cmd 


!  This  runtime  command  file  exercises  the  simulation  by  executing  load  changes.  In  this 
!  way  the  transient  response  of  a  source  buck  converter  can  be  analyzed.  The  file 
!  contains  a  study  that  allows  a  source  buck  converter  to  start  up.  The  final  values  of 
!  this  startup  are  used  as  ICs  for  the  transient  analysis.  In  this  way,  the  transient 
!  analysis  can  be  run  over  and  over  without  having  to  execute  a  startup  every  time. 

!  Written  by:  Michael  Badorf 
!  Last  modified:  02JUN97 

s  strplt  =  .t.  !"one  variable  per  x-axis  " 

s  calplt  =  .f. 


s  devplt  1  !"1  for  screen  output  " 

!  "5  for  postscript  " 

!  "true  rotates  plot  90  deg  " 

!  "x-axis  plot  units  " 

!  "false  suppresses  data  write  " 

!"each  time  SCHEDULE  occurs  " 

s  nrwitg  =  .t.  !  "true  enables  accumulation  of  data  " 

!  "after  a  CONTIN 

sdpnplt=.f.  !"no  date  on  plots  " 

!  "determines  which  variables  to  plot" 
prepare  t,iLl,ioutl,Vc,dutyl 


s  ppoplt  =  .£ 
s  xinspl  =  6 
s  weditg  =  .f. 


!  "startup  of  source  buck  converter" 
proced  sourcestartup 


.'"source  buck  startup  parameters" 
s  Rout  =  96.8  !"10%  load” 

!  "simulation  parameters" 
s  ialg  =  5 
s  maxt  =  1 ,0e-6 
s  cint  =  1.0e-3 
s  tstop  =  3.0 
start 


s  cint  =  1.0e-6 
s  tstop  =  4.0 
contin 


SAVE  /FILE-sourcestartup'  /BINARY  /STATES  /ICS  /EVENTS 
end 
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!  "single  source  buck  converter  transient  operation" 
proced  trans 

!  "initial  conditions  restoration  source  converter  startup" 
RESTORE  /FILE-sourcestartup'  /EVENTS 

["reinitializes  states  to  new  ICs" 

REINIT 

["source  buck  parameters" 

s  Rout  =  96.8  !"10%  load" 

!  "simulation  parameters" 
s  ialg  =  5 
s  maxt  =  1  .Oe-7 
s  cint  =  1.0e-6 
s  tstop  =  0.04 
start 

!"10%  ->  100%  source  buck  converter  load  change" 
s  Rout  =  9.68 
s  tstop  =  0.08 
contin 

!"100%  ->  10%  source  buck  converter  load  change" 
s  Rout  =  96.8 
s  tstop  =0.12 
contin 


splt=l 


plot  /color=0  /XTAG-  (sec)' 

& 

dutyl 

& 

ioutl  /TAG=,(amps)' 

& 

Vc  /TAG-(volts)’ 
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B.  LOAD  BUCK  CONVERTER  GAIN  SELECTION  CODE 

The  computer  code  for  this  section  is  comprised  of  the  main  ACSL  program  and 
the  command  file  used  to  generate  transient  load  studies.  The  buck  macro  is  provided 
above. 


1.  Main  Program 


!  loadbuck.csl 
! 

!  This  program  simulates  the  operation  of  a  3  kW  load  buck  converter  unit  operating  at 
!  20  kH2  assuming  ideal  components  and  no  line  losses.  The  converter  has  been  designed 
!  to  operate  as  follows: 

•  300  Vdc  input  voltage 

•  208  Vdc  output  voltage 

!  continuous  ind.  current  from  10%  ->  100%  load 

! 1116  controller  for  this  simulation  makes  use  of  the  digital  control  algorithm  derived 
!  in  Chapter  III.  The  file,  loadbuck.cmd,  contains  the  runtime  commands  utilized  to 
!  exercise  this  model. 

! 

!  Written  by:  Michael  Badorf 
!  Last  modified:  02JUN97 

IWM*****************,*****************,,*,**,**,,*,**,,*,*,,,,,,^^,,,,^,^^^^^^^^ 
INCLUDE  'buck.mac' 

PROGRAM  loadbuck 

INITIAL 


MAXTERVAL  maxt  =  1 ,0e-6 
MINTERVAL  mint  =  1.0e-8 


CINTERVAL  cint  =  1.0e-6 
ALGORITHM  ialg  =  5 


NSTEPS  nstp  =  1 
CONSTANT  tstop  =  1.0 

LOGICAL  reinitf 
CONSTANT  reinitf  =  .false. 


.'"maximum  integration  step  size" 

!  "minimum  integration  step  size" 
!"for  variable  step  algorithms  " 

!"data  communication  interval" 

!  "integration  algorithm" 

!"4~R.K.  2nd,  5~R.K.4th" 


.'"stop  point  for  integration" 
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|***************************************************************************** 
!  Initialization  of  digital  controller  values:  if  reinitializing  from  previous 
!  run,  then  this  section  is  skipped  by  setting  reinitf  =  .true,  in  the 
!  command  file. 

i ************************  ***************************************************** 

IF  (reinitf)  GO  TO  skipin 
Vthouse_old  =  0.0 
Verrold  =  0.0 
Vint_old  =  0.0 
count  =  0.0 
skipin.  .CONTINUE 


END  !  "of  initial” 
DYNAMIC 


TERMT  (t  .GE.  (tstop-0.5*cint)) 

DERIVATIVE 

Vincap3  =  300.0  !  "define  a  hard  input  voltage  of  300  Vdc" 

I* ***************************  ************************************************* 

!  load  buck  converter  unit  3 

I  ************* ***** ***** ** ********** ******* ***** ***** ******* **** * ************* 


!  "implement  load  buck  converter  unit  3  macro 

buck(3,Vincap3,iout3,isw3,iL3JVc3, count, "Lbuck3=1300.0e-6","Vdidrop3=2.0">& 
"Vswdrop3=2.0") 


CONSTANT  Cbuck3  =  400.0e-6 
CONSTANT  Vcic3  =0.0 

pVc3  =  (iL3  -  iout3)/Cbuck3 
Vc3  =  INTEG(pVc3,  Vcic3) 

CONSTANT  Rout3  =  144.0 


!  "output  capacitance  of  load  buck  " 

!  "output  voltage  IC  " 

!  "output  voltage  state  " 

!  "output  voltage  " 

!"load  buck  converter  startup  load  " 
!"10%  load  value  " 


iout3  =  Vc3/Rout3  !  "load  current" 

|*  ************************************************* 

!  end  load  buck  converter  unit  3 

j **************************************************** ************************* 


i  ***************************************************************************** 

!  load  buck  converter  unit  3  digital  controller  overcurrent  protection 
!************************** *************************************************** 


CONSTANT  iLlimit3  =  50.0 
CONSTANT  dutylow3  =  0.0 
CONSTANT  minduty3  =  0.0 
CONSTANT  maxduty3  =  475.0 
CONSTANT  scale3  =500.0 


!  "pulse-by-pulse  inductor  current  limit 
!  "lower  admissible  duty  cycle  limit 
!"lower  admissible  duty  cycle  limit 
!  "upper  admissible  duty  cycle  limit 
!"duty  between  0-500  counts 


t» 


tt 
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!  'if  load  buck  converter  operating,  implement  overcurrent  protection" 

IF  (iL3  .GT.  iLlimit3)  THEN 

duty3  =  dutylow3 !  "overcurrent" 

ELSE 

duty3  =  BOUND(minduty3,  maxduty3,  scale3*dutyML3) 

ENDIF  !"end  if  (iL3  .GT.  iLlimit3)" 

duty3p  =  duty3/scale3  !"scales  load  buck  duty  cycle  to  0  ->  1  range" 

!  end  load  buck  converter  unit  3  digital  controller  overcurrent  protection 

l*****************************************************^^^^*^^^^^** 
END  !  "of  derivative" 


!  implement  load  buck  converter  digital  controller 
DISCRETE  COUNTER 


!  "interval  statement  causes  controller  to  sample  every  tcount  seconds  " 

•'  tcount  large  if  load  buck  is  to  be  idle;  when  putting  on  line  " 

! "  reset  tcount  to  1 ,0e-7,  the  clock  speed  of  the  NSWC  digital  " 

!"  controller  ->  perform  reset  in  command  file 
INTERVAL  tcount  =  1 ,0e-7 

count  =  count  +  1  !  "dock  ->  advances  " 

!  "every  1.0e-7  seconds" 

CONSTANT  VthouserefB  =  208.0  !  "top  of  house  curve 

!  "reference  voltage  with  startup  "time  constant"  of  0.06  sec  " 

Vthouse3  =  Vthouse_old  +  700*tcount 

!  "memory  for  old  reference  voltage  value  -  used  on  startup  " 

Vthouseold  =  Vthouse3 


!  when  clock  counts  to  250  (—  25  usee)  update  duty  cycle  for " 
!  "application  to  upcoming  switching  interval  " 

IF  (count  .EQ.  250)  THEN 


!  "limits  reference  voltage  to  top  of  house  curve 
IF  (Vthouse3  .GE.  Vthouseref3)  THEN 
Vthouse3  =  Vthouseref3 
ENDIF  !"end  if  (Vthouse3  .GE.  Vthouseref3) 


CONSTANT  houseslope3  =  3.0 
CONSTANT  hi3  =  0.0105 
CONSTANT  hv3  =  0.000869 
CONSTANT  hn3  =  1.733 


!  "slope  of  house  curve 
!  "current  fdbk  gain 
!  "voltage  fdbk  gain 
.'"integral  voltage  fdbk  gain 
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!  "voltage  state  difference" 

Vprop  =  Vc3  -  Vthouse3 

! "integral  voltage  state  difference" 

Vint  =  25.0e-6*(Vprop  +  Verr_old)  +  Vint_old 

!  "control  algorithm  for  load  buck  converter  controller" 

dutysamp  =  -hi3*(iL3-iout3)-hv3*Vprop-hn3  Wint+Vthouse3/Vincap3 

Verr_old  =  Vprop 
Vint_old  =  Vint 

CALL  LOGD(.false.) 

ENDIF  !"end  if  (count  .EQ.  250)" 

!  "apply  duty  cycle  to  upcoming  switching  interval" 

IF  (count  .EQ.  500)  THEN 

CALL  LOGD(.false.) 

dutyML3  =  dutysamp 

count  =  0.0  ! "reset  clock  for  next  50  usee  cycle" 

CALL  LOGD(.false.) 

ENDIF  !"end  (count  .EQ.  500) 


END  !  "end  discrete 

END  !  "of  dynamic" 

END  !  "of  program" 
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2. 


Command  File 


loadbuck.cmd 


•  This  runtime  command  file  exercises  the  simulation  by  executing  load  changes.  In  this 
!  way  the  transient  response  of  a  load  buck  converter  can  be  analyzed.  The  file 
!  contains  a  study  that  allows  a  load  buck  converter  to  start  up.  The  final  values  of 
!  this  startup  are  used  as  ICs  for  the  transient  analysis.  In  this  way,  the  transient 
!  analysis  can  be  run  over  and  over  without  having  to  execute  a  startup  every  time. 

!  Written  by:  Michael  Badorf 
!  Last  modified:  02JUN97 

s  strplt  =  .t.  !  "one  variable  per  x-axis" 

s  calplt  = .{. 

s  devplt  =1  !  "1  for  screen  output" 

!  "5  for  postscript  " 


s  ppoplt  =  .f. 
sxinspl  =  6 
s  weditg  =  .f. 


!  "true  rotates  plot  90  deg” 

!  "x-axis  plot  units" 

!  "false  suppresses  data  write" 

!  "each  time  SCHEDULE  occurs" 


s  nrwitg  =  .t. 


!  "true  enables  accumulation  of  data" 
!  "after  a  CONTIN" 


s  dpnplt  =.f. 


!"no  date  on  plots" 


!  "determines  which  variables  to  plot" 
prepare  t,iL3,iout3,Vc3,duty3p 


!  "startup  of  load  buck  converter" 
proced  loadstartup 


!"load  buck  startup  parameters" 
s  Rout3  =  144.0  !"10%  load" 

!  "simulation  parameters" 
s  ialg  =  5 
s  maxt  =  1 ,0e-6 
s  cint  =  1.0e-6 
s  tstop  =  0.4 
start 


SAVE  /FILE-loadstartup’  count,  Vthouse_old,VeiT_old,Vint_old  & 
/BINARY /STATES /ICS /EVENTS 


end 
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["single  load  buck  converter  transient  operation" 
proced  trans 

["initial  conditions  restoration  from  load  buck  converter  startup" 
RESTORE  /FILE-loadstartup*  /EVENTS 

["reinitializes  state  to  new  ICs  and  skips  digital  controller" 

!  "initializations  in  INITIAL  section  of  main  program  " 

REINIT 
s  reinitf  =  .true. 

["load  buck  parameters" 

s  Rout3  =  144.0  ["10%  load" 

["simulation  parameters" 
s  ialg  =  5 
s  maxt  =  1.0e-7 
scint=  1.0e-6 
s  tstop  =  0.03 
start 

["10%  ->  100%  load  buck  converter  load  change" 
s  Rout3  =  14.4 
s  tstop  =  0.05 
contin 

["100%  ->  10%  load  buck  converter  load  change" 
s  Rout3  =  144.0 
s  tstop  =  0.08 
contin 


plot  /color=0  /XTAG=’(sec)’ 

& 

duty3p 

& 

ioutl  /TAG-Camps)’ 

& 

Vc  /TAG=,(volts)' 

end 
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C.  MODE  7  DETAILED  SIMULATION  CODE 

The  computer  code  for  this  section  is  comprised  of  the  main  ACSL  program,  three 
of  the  four  supporting  macro  files  and  the  command  file  used  to  generate  transient  load 
studies.  The  fourth  macro  file  is  the  analog  control  macro  provided  at  it  was  the  same  for 
all  simulation  cases. 


1.  Main  Program 


*****************t*************  *********************************^++^1^^^^ 

mode7.csl 


This  program  simulates  the  operation  of  parallel  source  buck  converter  units  (capable) 
of  supplying  1 8  k  W)  supplying  a  load  buck  converter  and  a  hard-switched  inverter 
The  buck  converter  models  assumes  equivalent  series  resistances  for  the 
capacitors  and  inductors.  The  inverter  operates  at  6  kHz.  The  system  has  been 
designed  to  operate  as  follows: 

Source  Buck  Converter:  400  Vdc  input  voltage  w/  360  Hz  ripple 
300  Vdc  nominal  output  voltage 
continuous  ind.  current  from  10%  ->  100%  load 
20  kHz  switching  frequency 

Load  Buck  Converter:  300  Vdc  nominal  input  voltage 

208  Vdc  nominal  output  voltage 
continuous  ind.  current  from  10%  ->  100%  load 
20  kHz  switching  frequency 


Inverter:  300  Vdc  nominal  input  voltage 

210  Vac(peak)  3-phase  output  voltage 

24  A(peak)  maximum  3-phase  current  output 

6  kHz  switching  frequency  (larger  value  requires  smaller  time  step) 


The  buck  converter  controllers  for  this  simulation  makes  use  of  the  analog  control 
algorithm  derived  in  Chapter  III.  The  inverter  is  controlled  by  stationary  reference 
frame  current  feedback  as  described  in  reference  [8],  The  file,  mode7.cmd,  contains 
the  runtime  commands  utilized  to  exercise  this  model. 


!  Written  by:  Michael  Badorf 
!  Last  modified:  02JUN97 

INCLUDE  'buck.mac' 
INCLUDE  'sourcefilter.mac' 
INCLUDE  'loadfilter.mac' 
INCLUDE  'analogcontrol.mac' 
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PROGRAM  operational  mode  7 
INITIAL 


MAXTERVAL  maxt  =  1.0e-6 
MINTERVAL  mint  =  1.0e-8 

CINTERVAL  cint  =  1.0e-6 
ALGORITHM  ialg  =  5 

NSTEPS  nstp=  1 
CONSTANT  tstop=  1.0 


!  "maximum  integration  step  size 


!"data  communication  interval 
!  "integration  algorithm 
!"4  -  R.K.  2nd,  5  -  R.K.  4th 

!"stop  point  for  integration" 


I***************************************************************************** 

!  Electrical  constants 

pi  =3.14159265 
twopi  =  2.0*pi 

pi2o3  =  2.0*3.14159265/3.0  !"  120  degrees  in  radians 

we  =60.0*2.0*pi  !  "60  Hz  converted  to  377  rad/sec  " 


!  Switch  and  diode  losses 

J ************************************** ********************* ******** ********** 

CONSTANT  Vdidrop  =  2.0  !  "voltage  drop  across  diodes  " 

CONSTANT  Vswdrop  =  2.0  !  "voltage  drop  across  IGBTs  " 

!  Buck  Converter  and  PWM  settings 

\ *********************************************%*********** ************ %%***%%* 
CONSTANT  swfreql  =  20000.0  {"switching  frequency  for  buck  1" 

Tbuckl  =  1.0/swfreql  {"switching  period  for  buck  1  " 


CONSTANT  swfreq2  -  20000.0  {"switching  frequency  for  buck  2" 
Tbuck2  =  1 .0/swfreq2  { "switching  period  for  buck  2  " 


CONSTANT  Vpk=  10.0 


i"pk  value  of  ramp  mod  wvfrm" 


CONSTANT  frnv  =  6000.0 
Tinv  =  1.0/finv 
Tinv_over2  =  Tinv/2.0 
slope  =  4.0*Vpk*fInv 


{"switching  frequency  for  inverter  ” 

!  "switching  period  of  inverter 
{"half  the  inverter  switching  period  " 
{"reference  triangle  slope  for  inverter" 


|******3tt*3(C*j|5******j|C***j|C***3|C3(c**3fCj|C*^**  +  **3jC5jc***3ic*)|C3|C*^**5|c*3jC3jt**jjt****3Jc****jJe^^*s|c*j|c*** 

{  Logical  switches  used  to  turn  load  buck  converter  on  and  off 

!  .  true  =  unit  on 
!  false  =  unit  off 

|*  ********  ******************  ************************* ****3^*********** ******** 

LOGICAL  buckl300 
CONSTANT  buckl  300  =  .false. 


119 


LOGICAL  reinitf 
CONSTANT  reinitf  =  .false. 

IF  (reinitf)  GO  TO  skipin 

I*************************************************]!******#******************** 
!  Initialization  of  digital  controller  values:  if  reinitializing  from  previous 
!  run,  then  this  section  is  skipped  by  setting  reinitf  =  .true,  in  the 
!  command  file. 

;****  *********  fc*******************************************************,,,,,..,,^.^* 

Vthouse3  =  0.0 
Verr_old  =  0.0 
Vintold  =  0.0 
count  =  0.0 

!  Initialization  of  inverter  values:  if  reinitializing  from  previous 
!  run,  then  this  section  is  skipped  by  setting  reinitf  =  .true,  in  the 
!  command  file. 

l********************************************^**^**************^^^^^^, 

Vag=0.0 

Vbg=0.0 

Vcg=0.0 

11  =  0.0 

12  =  0.0 

13  =  0.0 

skipin..CONTINUE 


END  !  "of  initial" 

DYNAMIC 

TERMT  (t  .GE.  (tstop  -  0.5*cint)) 
DERIVATIVE 


!  develops  6-pulse  voltage  input  to  source  buck  converters 
theta  =  we*t 

theta_s  =  mod(theta,  twopi)  ! "statement  causes  thetajs  to  vary  from  0  ->  2*pi  " 

!"then  reset  back  to  0  for  next  cycle  " 

CONSTANT  Vr  ave  =  208.0  !' "average  value  of  6  pulse  waveform  " 

Vrms  =  pi*Vr_ave/(3.0*SQRT(6.0))  !"rms  value  of  ripple  " 


!  "generates  6-pulse  waveform  " 

PROCEDURAL(V  r=theta_s,  Vrms,we) 

IF  ((theta_s  .GE.  O.O)  .AND.  (theta_s  .LE.  pi/3.0))  THEN 
Vr  =  -SQRT(6.0)*Vrms*cos(we*t  +  5.0*pi/6.0) 

ELSE  IF  ((theta_s  .GE.  pi/3.0)  .AND.  (theta_s  .LE.  2.0*pi/3.0))  THEN 
Vr  =  SQRT(6.0)*Vrms*cos(we*t  -  pi/2.0) 

ELSE  IF  ((theta_s  .GE.  2.0*pi/3.0)  .AND.  (theta_s  .LE.  pi))  THEN 
Vr  =  -SQRT(6.0)*Vrms*cos(we*t  +  pi/6.0) 

ELSE  IF  ((theta_s  .GE.  pi)  .AND.  (theta_s  .LE.  4.0*pi/3.0))  THEN 
Vr  =  SQRT(6.0)*Vrms*cos(we*t  +  5.0*pi/6.0) 

ELSE  IF  ((theta_s.GE.4.0*pi/3 .0). AND.(theta_s  .LE.  5.0*pi/3.0))  THEN 
Vr  =  -SQRT(6.0)*Vrms*cos(we*t  -  pi/2.0) 

ELSE  IF  ((theta_s  .GE.  5.0*pi/3.0)  .AND.  (theta_s  .LE.  2.0*pi))  THEN 
Vr  =  SQRT(6.0)*Vrms*cos(we*t  +  pi/6.0) 

ENDIF  !"end  multiple  if-then-else" 

END  !  "of  procedural" 

Vin  =  2*Vr  !  "voltage  at  source  buck  converter  input  filters" 

I***************************************************************************** 

!  end  6-pulse  voltage  input  to  source  buck  converters 

|**3fc****9|e*3jc^********4c*^%*)|e^4e*^*3i(4nie^*3fe**^]tc9ic}|c**4;$ic*>tc3|(}|c***^**«4:3|e*t**************** 

!  source  buck  converter  unit  1 

i ***************************************************************************** 

!  "input  LC  filter  for  source  buck  converter  unit  1” 

sourcefilter(  1 ,  Vin,  Vincap  1  ,iLfilt  1  ,isw  1  ,"Rlesr_in  1  =0.2",  "Lfilt 1  =425 .0e-6",& 

"Rcesr_inl  =2046.6", "Cin  1  =2000.0 1  e-6","Vcapic  1  =400.0") 

tpwml  =  MOD(t,Tbuckl)  !"time  elapsed  in  switching  period  for  buck  " 

Vramp  1  =  Vpk*tpwm  1/Tbuck  1  !  "establish  ramp  modulating  waveform  for  buck 

["implement  source  buck  converter  unit  1  macro” 

buck(  1 ,  Vincap  1  ,iout  1  ,isw  1  ,iL  1  ,Vout,  Vramp  1  ,”Lbuck 1 =760.0e-6",& 

"  V  didropl  =2.0","  V  swdrop  1  =2 .0") 


["implement  source  buck  converter  unit  1  multiloop  control 
analog_control(  1 ,  Vincap  1 ,  Vout,iout  1  ,iL  1  ,"RCstart  1  =0. 1", "hi  1  =0.0 1 5 
"hv  1  =0.0 1 7", "hn  1 =26.09") 

!  end  source  buck  converter  unit  1 

!  source  buck  converter  unit  2 

|5k*******3}e***^**)k**5t:*****^******5jC*5|t  +  **5»:3je5tC****3|C**  +  *  +  ***^3({3(t^^^^^*^}jl:3)C5t-^:4c  +  ***^^*  +  ^ 

["input  LC  filter  for  source  buck  converter  unit  2” 
sourcefilter(2,Vin,Vincap2,iLfilt2,isw2,"Rlesr_in2=0.2","Lfilt2=425.0e-6",& 
”Rcesr_in2=2046.6","Cin2=2000.0 1  e-6”,"  Vcapic2=400.0”) 
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tpwm2  M0D(t,Tbuck2)  !”time  elapsed  in  switching  period  for  buck  " 

Vramp2  =  Vpk*tpwm2/Tbuck2  !  "establish  ramp  modulating  waveform  for  buck  " 

!  "implement  source  buck  converter  unit  2  macro” 

buck(2,Vincap2,iout2,isw2,iL2,Vout,Vramp2,"Lbuck2=875.0e-6"  & 
"Vdidrop2=2.0","Vswdrop2=2.0") 

■  implement  source  buck  converter  unit  2  multiloop  control  macro” 

analog_control(2,Vincap2,Vout,iout2,iL2,"RCstart2=0.r’,"hi2=0.015",& 

"hv2=0.0 1 7'V'hn2=26.09") 

!  end  source  buck  converter  unit  2 

l*****************************************************^^^^^^^^ 


CONSTANT  Cbuck  =  810.0e-6 
CONSTANT  Rout  =  96.8 
CONSTANT  Rcesr  =  0.02 
CONSTANT  Vcic  =  0.0 

pVc  =  (iLl+iL2-iout)/Cbuck 
Vc  =  INTEG(pVc,  Vcic) 

Vout  =  Vc  +  (iL  1  +iL2-iout)*Rcesr 


ioutl  =  iLl  -  (Cbuck/2)*pVc 
iout2  =  iL2  -  (Cbuck/2)*pVc 


!  "output  capcitance 

.'"source  buck  converter  startup  load 

["capacitor  esr 

["output  voltage  IC 

!  "output  capacitor  voltage  state  " 

!  "voltage  across  output  capacitor  " 

["output  voltage  " 

["total  iout 

["ioutofbuck  1  " 

!"iout  of  buck  2  " 


iout  -  (Vc+(iLl+iL2)*Rcesr+(iLfilt3+iinv)*Rout)/(Rout+Rcesr) 


ft 


M 

M 


M 


load  buck  converter  unit  3 


["input  LC  filter  for  load  buck  converter  unit  3” 

loadfilter(3,Vout,Vincap3,iLfilt3,isw3,"Rlesr_in3=0.2","Lfilt3=425.0e-6",& 

"Rcesr_in3=0.01","Cin3=460.0e-6","Vcapic3=0.0") 

["implement  source  buck  converter  unit  3  macro” 

buck(3,Vincap3,iout3,isw3,iL3,Vout3, count, "Lbuck3=1300.0e-6","VdidroD3=2  0"  & 

"Vswdrop3=2.0")  ’  ’ 


CONSTANT  Cbuck3  =  400.0e-6 
CONSTANT  Vcic3  =  0.0 
CONSTANT  Rout3  =  144.0 
CONSTANT  Rcesr_out3  =  0.01 


pVc3  =  (iL3  -  iout3)/Cbuck3 
Vc3  =  INTEG(pVc3,  Vcic3) 


.'"output  capcitance 

["output  voltage  IC 

["source  buck  converter  startup  load 

["load  buck  converter  output 

["capacitor  esr 

["capacitor  output  voltage  state 
["capacitor  output  voltage 


Vout3  -  Vc3  +  (iL3-iout3)*Rcesr_out3  [’’output  voltage 

iout3  =  (Vc3+iL3*Rcesr_out3)/(Rout3+Rcesr_out3)  •  ["total  iout 


it 

It 


ft 


tf 

It 


tt 
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[***  ************  ************************************************************** 
!  end  load  buck  converter  unit  3 

I*** *********** *************************************************************** 


!  load  buck  converter  unit  3  digital  controller  overcurrent  protection 

i************************************************************** *************** 


CONSTANT  iLlimitS  -  50.0 
CONSTANT  dutylow3  =  0.0 
CONSTANT  minduty3  -  0.0 
CONSTANT  maxduty3  =  475.0 
CONSTANT  scale3  =500.0 


!  "pulse-by-pulse  inductor  current  limit " 
!  "lower  admissible  duty  cycle  " 

!  "lower  admissible  duty  cycle  limit  " 

!  "upper  admissible  duty  cycle  limit  " 
!"duty  between  0-500  counts" 


IF  (buck 13  00)  THEN 


!  "if  load  buck  converter  operating,  implement  overcurrent  protection" 
IF  (iL3  .GT.  iLlimit3)  THEN 


duty3  =  dutylow3!  "overcurrent" 


ELSE 


duty3  =  BOUND(minduty3,  maxduty3,  scale3*dutyML3) 
ENDIF  !"end  if  (iL3  .GT.  iLlimitS)" 


ELSE 

! "if  load  buck  converter  is  idle,  duty  cycle  is  zero" 
duty3  =  0.0 
Vout3  =  0.0 
iout3  =  0.0 

ENDIF  !  "end  if  (buck  1 300)" 

duty3p  =  duty3/scale3  !  "scales  load  buck  duty  cycle  to  0  ->  1  range" 

l***************************************************************************** 

!  end  load  buck  converter  unit  3  digital  controller  overcurrent  protection 

l***************************************************************************** 


!  inverter  model 


CONSTANT  Lf 

=  10.0e-3 

!  "inverter  phase  load  inductance 

CONSTANT  R4 

=  5.0 

!  "inverter  phase  load  resistance 

CONSTANT  K 

=  0.0005 

!"PWM  gain 

CONSTANT  K1 

=  1260.0 

!  "proportional  gain 

CONSTANT  K2 

=  3750.0 

!  "integral  gain 

CONSTANTA 

=  24.0 

!  "desired  current  amplitude 
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!  "sine  -  triangle  modulation" 
tt  =  mod(t,Tinv) 

IF(tt  .LT.  Tinv_over2)  THEN 

Vtri  =  Vpk  -  tt*slope 

ELSE 

Vtri  =  -Vpk  +  slope*(tt  -  Tinv_over2) 

ENDIF  !"end  if  (tt  .LT.  Tinv_over2)" 

!  "calculate  synchronous  reference  frame  angle" 
thetae  =  integ(we,0.0) 

!  "calculate  desired  currents" 

Idesa  =  A*cos(thetae)  !  "desired  phase  A  current  " 

Idesb  =  A*cos(thetae  -  pi2o3)  .'"desired  phase  B  current  " 

Idesc  =  A*cos(thetae  +  pi2o3)  ("desired  phase  C  current  " 

("transform  desired  currents  to  stationary  qdO  quantities" 

Iqsdes  =  (2.0/3.0)*(Idesa  -  0.5*Idesb  -  0.5*Idesc) 

Idsdes  =  (sqrt(3.0)/3.0)*(  -Idesb  +  Idesc) 

("transform  desired  currents  to  synchronous  qdO  quantities" 

Iqsedes  =  (2.0/3.0)*(Idesa*cos(thetae)  +  Idesb*cos(thetae  -  pi2o3)  & 

+  Idesc*cos(thetae  +  pi2o3)) 

Idsedes  =  (2.0/3.0)*(Idesa*sin(thetae)  +  Idesb*sin(thetae  -  pi2o3)  & 

+  Idesc*sin(thetae  +  pi2o3)) 

!  transform  actual  currents  to  synchronous  qdO  quantities" 

Iqse  =  (2.0/3.0)*(ILa*cos(thetae)  +  ILb*cos(thetae  -  pi2o3)  & 
+ELc*cos(thetae  +  pi2o3)) 

Idse  =  (2.0/3 .0)*  (ILa*sin(thetae)  +  ILb*sin(thetae  -  pi2o3)  & 
+ILc*sin(thetae  +  pi2o3)) 

("determine  error  currents  in  stationary  reference  frame” 

Iqserr  =  Iqsdes  -  iqss 
Idserr  =  Idsdes  -  idss 

("generate  commanded  currents  (PI  Controller)" 

Iqscom  =  (K*K1  *Iqserr  +  K*K2*integ(Iqserr,0.0)) 

Idscom  =  (K*Kl*Idserr  +  K*K2*integ(Idserr,0.0)) 

("transform  commanded  quantaties  to  ABC  frame" 

Iacont  =  Iqscom 

Ibcont  =  -0.5*Iqscom  -  (sqrt(3.0)/2.0)*Idscom 
Iccont  =  -0.5*Iqscom  +  (sqrt(3.0)/2.0)*Idscom 
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{’’switching  control  voltages" 

Vacont  =  Iacont*R4 
Vbcont  =  Ibcont*R4 
Vccont  =  Iccont*R4 

!  "schedule  inverter  switching" 

SCHEDULE  hia  .XN.  (Vtri  -  Vacont) 

SCHEDULE  loa  .XP.  (Vtri  -  Vacont) 

SCHEDULE  hib  .XN.  (Vtri  -  Vbcont) 

SCHEDULE  lob  .XP.  (Vtri  -  Vbcont) 

SCHEDULE  hie  .XN.  (Vtri  -  Vccont) 

SCHEDULE  loc  .XP.  (Vtri  -  Vccont) 

{"calculate  line  voltages" 

Vas  =  (2.0/3 .0)*Vag  -  (1.0/3.0)*(Vbg  +  Vcg) 

Vbs  =  (2.0/3 .0)*Vbg  -  (1.0/3.0)*(Vag  +  Vcg) 

Vcs  =  (2.0/3.0)*Vcg  -  (1.0/3.0)*(Vag  +  Vbg) 

{"transform  to  stationary  reference  frame" 

Vqss  -  (2.0/3. 0)* (Vas  -  0.5*Vbs  -  0.5*Vcs) 

Vdss  -  (sqrt(3.0)/3.0)*(  -Vbs  +  Vcs) 

{"find  current  derivatives  in  stationary  frame" 
piqss  -  (Vqss  -  R4*iqss)/Lf 
pidss  -  (Vdss  -  R4*idss)/Lf 

{"calculate  stationary  currents" 

CONSTANT  iqssic  =  0.0  {"initial  q  axis  current  " 

CONSTANT  idssic  =  0.0  {"initial  d  axis  current  " 

iqss  —  integ(piqss5iqssic) 
idss  =  integ(pidss?idssic) 

!  "transform  back  to  ABC  quantaties" 
iLa  =  iqss 

iLb  =  -0.5*iqss  -  (sqrt(3.0)/2.0)*idss 
iLc  =  -0.5*iqss  +  (sqrt(3.0)/2.0)*idss 

iinv  =  il  +  i2  +  i3 

END  { "of  derivative" 

\***************************************************************************** 

{  implement  load  buck  converter  digital  controller 
|  **  *****  *  *  *  *  ***************************  ***  **  **  ******  *  ***  *  *****^*  +  ****^^  +  +  *  +  *^  * 

DISCRETE  COUNTER 

{"interval  statement  causes  controller  to  sample  every  tcount  seconds  " 

!"  set  tcount  large  if  load  buck  is  to  be  idle;  when  putting  on  line" 

{"  reset  tcount  to  1  .Oe-7,  the  clock  speed  of  the  NSWC  digital  " 

{ "  controller  ->  perform  reset  in  command  file 
INTERVAL  tcount  =  10.0 


{"switch  A  to  HI  state 
{"switch  A  to  LO  state 
{"switch  B  to  HI  state 
{"switch  B  to  LO  state 
{"switch  C  to  HI  state 
{"switch  C  to  LO  state 
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!  "execute  if  block  when  load  buck  converter  is  to  be  put  on  line  " 

! "  the  buck 1 3 00  condition  tells  the  program  that  the  load  buck 
! "  is  to  be  operating;  the  tcount  condition  ensures  that  tcount " 

!"  has  been  reset  since  reset  only  takes  effect  on  following  " 

•"  program  cycle  ->  see  ACSL  Reference  Manual  for  details  " 

IF  ((buckl300)  .AND.  (tcount  .EQ.  1.0e-7))  THEN 

count  =  count  +  1  !  "clock  ->  advances  " 

!  "every  1.0e-7  seconds" 

CONSTANT  Vthouserefi  =  208.0  !"top  of  house  curve  " 

!"reference  voltage  with  startup  "time  constant"  of  0.06  sec  " 

Vthouse3  =  Vthouseold  +  700*tcount 

["memory  for  old  reference  voltage  value  -  used  on  startup  " 

Vthouseold  =  Vthouse3 

!  when  clock  counts  to  250  (=  25  usee)  update  duty  cycle  for  " 

!  "application  to  upcoming  switching  interval  " 

IF  (count  .EQ.  250)  THEN 

["limits  reference  voltage  to  top  of  house  curve  " 

IF  (Vthouse3  .GE.  Vthouseref3)  THEN 
Vthouse3  =  Vthouseref3 
ENDIF  ["end  if(Vthouse3  .GE.  Vthouseref3) 

CONSTANT  houseslope3  =  3.0  ["slope  of  house  curve 

CONSTANT  hi3  =  0.0105  ["current  fdbk  gain  " 

CONSTANT  hv3  =  0.000869  ["voltage  fdbk  gain  " 

CONSTANT  hn3  =  1.733  ["integral  voltage  fdbk  gain" 

["voltage  state  difference" 

Vprop  =  Vout3  -  Vthouse3 

["integral  voltage  state  difference" 

Vint  =  25.0e-6*(Vprop  +  Verr_old)  +  Vint_old 

["control  algorithm  for  load  buck  converter  controller" 

dutysamp  =  -hi3*(iL3-iout3)-hv3*Vprop-hn3*Vint+Vthouse3/Vincap3 

Verr_old  =  Vprop 
Vintold  =  Vint 

CALL  LOGD(.false.) 

ENDIF  !  "end  if  (count  .EQ.  250)" 
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!  "apply  duty  cycle  to  upcoming  switching  interval  " 

IF  (count  .EQ.  500)  THEN 

CALL  LOGD(.false.) 

dutyML3  =  dutysamp 

count  =  0.0  ["reset  clock  for  next  50  usee  cycle" 

CALL  LOGD(.false.) 

ENDIF  !"end  (count  .EQ.  500) 

ENDIF  !”endif(buckl300)...  " 

END  !"end  discrete  " 

DISCRETE  hia 

Vag  =  Vout  -  Vswdrop 
il  =  iLa 

END  !  "of  discrete" 

DISCRETE  loa 

Vag  =  Vswdrop 

11  =0.0 

END  !  "of  discrete" 

DISCRETE  hib 

Vbg  =  Vout  -  Vswdrop 

12  =  iLb 

END  !  "of  discrete" 

DISCRETE  lob 

Vbg  =  Vswdrop 

12  =0.0 

END  !  "of  discrete" 

DISCRETE  hie 

Vcg  =  Vout  -  Vswdrop 

13  =  iLc 

END  !  "of  discrete" 

DISCRETE  loc 

Vcg  =  Vswdrop 
i3  =0.0 

END  !  "of  discrete" 

END  !  "of  dynamic" 

END  !  "of  program” 
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2. 


Buck  Converter  Macro 


!  buck.mac 
! 

!  This  macro  governs  the  conduction  status  of  a  buck  converter  switch  and  uses  this 
!  condition  to  develop  inductor  current  utilizing  the  inductor  equation  of  state. 

!  Switch  and  diode  losses  are  modeled  as  voltage  drops.  Capacitors  and  inductors  are 
!  modeled  with  equivalent  series  resistance  in  place.  The  output  inductor  current  is 
!  used  in  the  main  program  to  determine  output  capacitor  voltage.  Although  this  step 
!  could  be  done  in  the  macro,  the  current  arrangement  allows  modeling  more  complex  PEBB 
!  network  configurations. 

i 

!  Written  by:  Michael  Badorf 
!  Last  modified:  02JUN97 

1********************************************************^^***^^^*^********* 

MACRO  buck(z,vin,iout,isw,iL,Vo,rampwav,pLbuck,pVdidrop,pVswdrop) 


INITIAL 


CONSTANT  pLbuck 
CONSTANT  pVdidrop 
CONSTANT  pVswdrop 
CONSTANT  iLic&z  =  0.0 
CONSTANT  Rlesr&z  =  0.2 
CONSTANT  Rcesr&z  =  0.01 


"inductance  " 
"diode  losses  " 

"switch  losses  " 

"inductor  current  IC" 
"inductor  esr  " 
"capacitor  esr  " 


"——define  and  initialize  the  switch  conduction  status" 
LOGICAL  SWl&z 

SWl&z  =  .true.  !"T  si  conducts,  F  si  blocks" 


END  !”of  initial” 


" — determine  conduction  status  of  SW1 " 
PROCEDURAL(SW  1  &z,isw=rampwav,duty&z,iL) 

" — charge  inductor  when  duty  GT  reffamp" 
IF  (duty&z  .GT.  rampwav)  THEN 
SWl&z  =  .true. 
isw  =  iL 

ELSE 

SWl&z  =  . false, 
isw  =  0.0 

ENDIF 

END  !  "of  procedural" 
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" — Given  conduction  status  of  SW1 ,  and  whether  " 

"  contin  or  not,  determine  deriv  of  iL  " 

PROCEDURAL(piL&z=SWl&z,vin,Vo,Vdidrop&z,iL,iout,Vswdrop&z,Lbuck&z,& 

RJesr&z,Rcesr&z) 

" — if  SW1  gated,  assume  can  conduct  pos.  current" 

IF  (SWl&z)  THEN 

piL&z  =  (vin-Vswdrop&z-Vo)/Lbuck&z  -  & 
(Rlesr&z+Rcesr&z)*iL/Lbuck&z  +  & 
Rcesr&z*iout/Lbuck&z 

ELSE 

IF  (iL  .GT.  0.0)  THEN 

!"if  IGBT  not  gated  but  iL  greater  than  zero  " 

!"the  diode  conducts  and  the  inductor  discharges" 
piL&z  =  (-Vdidrop&z-Vo)/Lbuck&z  -  & 

(RIesr&z+Rcesr&z)*iL/Lbuck&z  +  & 

Rcesr&z*  iout/Lbuck&z 

ELSE 

!"if  IGBT  not  gated  and  iL  less  than  or  equal" 

!"to  zero,  the  diode  is  off  and  the  inductor  " 

!  "current  and  derivative  are  zero  " 

piL&z  =  0.0 

ENDIF  !"end  if  (iL  .GT.  0.0)" 

ENDIF  !"end  if  (SWl&z)" 

END  !  "of  procedural" 

!  "integrate  the  inductor  state  variable" 

iL  =  BOUND(0.0,  100.0,  LIMINT(piL&z,  iLic&z,  0.0,  100.0)) 

MACRO  END 


I 
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3. 


Source  Buck  Converter  Input  Filter  Macro 


ji«**iMf*******##**#iM****M******iM*************************************************** 
!  sourcefilter.mac 
! 

!  This  macro  governs  the  operation  of  the  source  buck  converter  input  LC  filters.  The 
!  model  takes  into  account  equivalent  series  resistance  of  the  capacitor  and  inductor  in 
!  the  filter.  Due  to  algebraic  loop  considerations  within  the  controller  and  buck 
!  macros,  an  equivalent  parallel  capacitance  and  resistance  was  found  for  the  series 
!  capacitance/esr  combination.  This  result  gives  the  seemingly  large  value  of  capacitor 
!  esr  shown  in  the  call  to  this  macro  in  mode7.csl. 

i 

!  Written  by:  Michael  Badorf 
!  Last  modified:  02JUN97 

MACRO  sourcefilter(z,vin,Vcap,iLfilt,isw,pRlesr_in,pLfilt,pRcesr_in,& 

pCin,pVcapic) 

INITIAL 

CONSTANT  pRlesrJn 
CONSTANT  pLfilt 
CONSTANT  pRcesrin 
CONSTANT  pCin 
CONSTANT  pVcapic 

" — assign  the  inductor  current  IC" 

CONSTANT  iLfiltic&z  =  0.0 

END  !  "of  initial" 

!  "buck  converter  input  filter  state  equations  " 
piLfilt&z  =  (vin-Rlesr_in&z*iLfilt-Vcap)/Lfilt&z 

!  "current  through  input  inductor  " 
iLfilt  =  INTEG(piLfilt&z,  iLfiltic&z) 

pVcap&z  =  (iLfilt  -  isw  -  Vcap/Rcesr_in&z)/Cin&z 

!  "input  voltage  at  buck  converter  IGBT  " 

Vcap  =  INTEG(pVcap&z,  Vcapic&z) 

MACRO  END 
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4. 


Load  Buck  Converter  Input  Filter  Macro 


| *** *********************************************************************** *********** 
!  loadfilter.mac 
! 

!  This  macro  governs  the  operation  of  the  load  buck  converter  input  LC  filter.  The 
!  model  takes  into  account  equivalent  series  resistance  of  the  capacitor  and  inductor  in 
!  the  filter.  In  contrast  to  the  source  filter,  no  algebraic  loops  were  formed  and  the 
!  series  capacitance/esr  could  be  modeled  directly. 

I 

!  Written  by:  Michael  Badorf 
!  Last  modified:  02JUN97 

I* ************************** ***************** ***************************************** 

MACRO  loadfilter(z,vin, V cap, iLfilt, isw,pRlesr_in,pLfilt,pRcesr_in,& 
pCin,pVcapic) 


INITIAL 

CONSTANT  pRlesr  in 
CONSTANT  pLfilt 
CONSTANT  pRcesrin 
CONSTANT  pCin 
CONSTANT  pVcapic 

" — assign  the  initial  inductor  current  and  cap  volt” 
CONSTANT  iLfiltic&z  -  0.0 


END  !"of  initial" 


!"buck  converter  input  filter  state  equations  " 
piLfilt&z  =  (vin-Rlesr_in&z*iLfilt-Vcap)/Lfilt&z 

!  "current  through  input  inductor 
iLfilt  =  INTEG(piLfilt&z,  iLfiltic&z) 

pV&z  =  (iLfilt  -  isw)/Cin&z 
V&z  -  INTEG(pV&z,  Vcapic&z) 

!  "input  voltage  at  buck  converter  IGBT  " 

Vcap  =  V&z  +  (iLfilt-isw)*Rcesr_in&z 

MACRO  END 
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5. 


Command  File 


*********************************^**^*****#******************>1^************^*^**#** 

mode7.cmd 


!  This  runtime  command  file  exercises  the  simulation  by  executing  load  changes.  In  this 
!  way  the  transient  response  of  the  PEBB  network  can  be  analyzed.  The  file  contains 
!  studies  that  allow  system  start  up.  The  final  values  of  the  startups  are  used  as  ICs 
!  for  the  transient  analysis.  In  this  way,  the  transient  analysis  can  be  run  over  and 
!  over  without  having  to  execute  a  startup  every  time. 

i 

!  Written  by:  Michael  Badorf 
!  Last  modified:  02JUN97 

l**********************************#*^*^^*****************^^*************^*^** 

s  strplt  =  .t.  !  "one  variable  per  x-axis" 

s  calplt  =  .f. 


s  devplt  =  5 


s  ppoplt  = .{. 
s  xinspl  =  6 
s  weditg  =  .f. 


!  "6  for  X-windows" 

!  "5  for  postscript" 

!  "true  rotates  plot  90  deg" 

!  "x-axis  plot  units" 

!  "false  suppresses  data  write" 

!  "each  time  SCHEDULE  occurs" 


s  nrwitg  =  .t. 


!  "true  enables  accumulation  of  data" 
!  "after  a  CONTIN" 


s  dpnplt  =  .f. 


!  "no  date" 


!  "determines  which  variables  to  plot" 

prepare  t,iout,ioutl,iout2,iout3,Vout,Vout3,iLa,Vas 

!  "startup  of  source  buck  converters" 
proced  startup 


!  "source  buck  startup  parameters" 
s  Rout  =  40.0 


Hoad  buck  parameters" 
s  Rout3  =  144.0 

sbuckl300  =  . false.  !"loadbuck  off’ 

s  tcount  =  10.0 

!  "inverter  parameters 
s  R4  =  5.0 

s  A  =  0.0  !  "inverter  off' 
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["simulation  parameters" 
s  ialg  =  5 
s  maxt  =  1  .Oe-6 
s  cint  =  1  .Oe-3 
s  tstop  =  0.5 
start 

s  cint  =  1  .Oe-6 
s  tstop  =  0.75 
contin 

SAVE  /FILE- startup'  count, Vthouse_old,Verr_old,Vint_old,& 
Vout,iout,  V  incap3 ,  Vout3  ,iout3 ,  Vag,  Vbg,  Vcg,  il,i2,i3,& 
/BINARY  /STATES  /ICS  /EVENTS 


!  "startup  of  load  buck  converter" 
proced  loadbuck 

RESTORE  /FILE-startup’  /EVENTS 
s  reinitf  =  .true. 

REINIT 

["source  buck  startup  parameters" 
s  Rout  =  40.0 

!"load  buck  parameters" 
s  Rout3  =  144.0 

s  buckl300  -  .true.  ! "loadbuck  on" 

stcount  =  1.0e-7 

!  "inverter  parameters 
s  R4  =  5.0 

s  A  =  0.0  ["inverter  off' 

!  "simulation  parameters" 
s  ialg  =  5 
smaxt  =  1.0e-6 
s  cint  =1. Oe-6 
s  tstop  =  0.4 
start 

SAVE  /FILE- loadbuck’  count, Vthouse_old,Verr_old,Vint_old,& 
Vout,iout,Vincap3,Vout3,iout3,Vag,Vbg,Vcg,il,i2,i3,& 
/BINARY /STATES /ICS /EVENTS 


end 
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.'"startup  of  inverter" 
proced  inverter 

RESTORE  /FILE- loadbuck'  /EVENTS 
s  reinitf  =  .true. 

REINIT 

!  "source  buck  startup  parameters" 
s  Rout  =  40.0 

Hoad  buck  parameters" 
s  Rout3  =  144.0 

s  buck  1300  =  .true.  !  "loadbuck  on" 

s  tcount  =  1.0e-7 

.'"inverter  parameters 
s  R4  =  5.0 

s  A  =  10.0  !"inverter  on" 


!  "simulation  parameters" 
s  ialg  =  5 
s  maxt  =  1.0e-6 
s  cint  =  1.0e-6 
s  tstop  =  0.4 
start 

SAVE  /FILE=’inverter'  count, Vthouse_old,Verr_old,Vint_old,& 
V  out,iout,  V  incap3,  V  out3,iout3 ,  Vag,  Vbg,  Vcg,i  1  ,i2,i3  ,& 
/BINARY  /STATES  /ICS  /EVENTS 


end 

.'"parallel  source  buck  converters  supplying  load  buck  and  inverter" 
proced  trans 

.'"initial  conditions  restoration  from  760  uH  source  " 

!  "buck  converter  /  load  buck  converter  startup  " 
RESTORE  /FILE-inverter'  /EVENTS 
REINIT 
s  reinitf  =  .true. 

!  "source  buck  parameters" 
s  Rout  =  96.8 

Hoad  buck  parameters" 
s  Rout3  =  144.0 
s  buck 1300  =  .true, 
s  tcount  =  1 ,0e-7 

!  "inverter  parameters" 
s  R4  =  5.0 
sA  =  10.0 
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!  "simulation  parameters" 
s  nrwitg  =  .true, 
s  ialg  =  5 
smaxt  =  1.0e-6 
s  cint  =  1.0e-6 
s  tstop  =  0.04 
start 

!"10%  ->  100%  load  buck  converter  load  change" 

!"3  kW  100%  power  reference  " 

s  Rout3  =  14.4 
s  tstop  =  0.08 
contin 

!"100%  ->  10%  load  buck  converter  load  change" 
s  Rout3  -  144.0 
s  tstop  =0.12 
contin 

!"40%  ->  100%  inverter  load  change" 
s  A  =  24.0 
s  tstop  =  0.17 
contin 

!"100%  ->  40%  inverter  load  change” 
s  A  =  10.0 
s  tstop  =  0.2 
contin 

s  plt=l 

plot  /color=0  /XTAG-(sec)’  & 

iout3  /TAG=’(amps)'  /lo=0.0  /hi  =16.0  & 

Vout3  /TAG-(volts)' /lo=l 96.0 /hi=2 18.0  & 
iout  /TAG- (amps)’  /lo=-5.0 /hi=30.0  & 

Vout  /TAG- (volts)'  /lo=300.0  /hi=308.0 


s  plt=2 

plot  /color=0  /XTAG=,(sec)'  & 

iout2 /TAG- (amps)’  /lo=-5.0  /hil6.0  & 

ioutl  /TAG- (amps)’  /lo=-5.0  /hi=16.0  & 

iLa  /TAG-  (amps)1  /lo=-24.0 /hi=24.0  & 

Vas  /TAG=’(volts)’/lo=-210.0/hi=210.0 


end 
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APPENDIX  D.  EASYTRAX  PCB  DATA 


A.  ANALOG  CONTROLLER  PCB  DATA 

This  section  is  divided  into  a  component  parts  listing,  an  EASYTRAX  PCB  netlist,  and  a  PCB 
component  overlay. 


1.  Component  List 


R411G 

AXIAL0.4 

AXIAL0.4 

R807Q 

C603 

AXIAL0.4 

RAD0.2 

R807H 

R603 

AXIAL0.4 

AX3AL0.4 

R606 

AX3AL0.4 

R603H 

AXIAL0.4 

R803 

AXIAL0.4 

R411 

AX3AL0.4 

R803G 

AXIAL0.4 

R607 

AXIAL0.4 

C107 

RADO.2 

C902 

RADOO.2 

D107 

DIODEO.4 

C904 

RB.3/.6 

D106 

DIODEO.4 

C905 

RB.3/.6 

R112G 

AXIAL0.4 

C907 

RB.2.4 

R114 

AXIAL0.4 

C901 

RB.2.4 

R108 

AXIAL0.4 

R112 

R110G 

C203 

AXIAL0.4 

RADO.2 

R109 

C210 

AXIAL0.4 

RADO.2 

R102B 

C212 

AXIAL0.4 

RADOO.2 

R113A 

C205 

AX3AL0.4 

RADO.2 

R113C 

R024J 

AXIAL0.4 

AXIAL0.4 

R106 

AXIAL0.4 

R205 

R313 

AXIAL0.4 

AXIAL0.4 

R210 

R102C 

AXIAL0.4 

R212 

AXIAL0.4 

AXIAL0.4 

R101 

C904D 

AXIAL0.4 

RADO.2 

R113B 

C305 

AXIAL0.4 

RADO.2 

R110 

C804 

AXIAL0.4 

RADO.2 

R102A 

C808 

AXIAL0.4 

RADO.2 
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C514 

RAD0.2 

C708 

RAD0.2 

C704 

RAD0.2 

C306 

RAD0.2 

cm 

RAD0.2 

C104 

RAD0.2 

C21 1 
RAD0.2 

C204 

RAD0.2 

C608 

RAD0.2 

C414 

RAD0.2 

C408 

RAD0.2 

C604 

RAD0.2 

D901A 

DIODEO.4 

D902A 

DIODEO.4 

D901B 

DIODEO.4 

D902B 

DIODEO.4 


R902 

AXIAL0.4 

R907 

AXIAL0.4 

D907 

RADO.l 

R313G 

AXIAL0.4 

R704 

AXIAL0.4 

R702 

AXIAL0.4 

R316B 

AXIAL0.4 

R316A 

AXIAL0.4 

R318 

AXIAL0.4 

R301 

AXIAL0.4 

R311 

AXIAL0.4 

C313 

RAD0.2 

C301 

RAD0.2 

C311 

RAD0.2 

C302 

RAD0.2 

C317 

RAD0.2 


C702 

RAD0.2 


D311 

DIODEO.4 

D301 

DIODEO.4 

D317 

DIODEO.4 

R501 

AXIAL0.4 

R509 

AXIAL0.4 

R504 

AXIAL0.4 

R807 

AXIAL0.4 

R505 

AXIAL0.4 

R502 

AXIAL0.4 

C807 

RAD0.2 

C501 

RAD0.2 

C803 

RAD0.2 

C316 

RAD0.2 

R603V 

VR1 


J25 

DB25RA/F 

C907G 

RAD0.2 

U1 

DIP  14 
LF347 

U2 

DIP  12 
LM324 

U3 

DIP  18 
UC3637 

U4 

DIP  14 
AD534 

U5 

DIP  14 
7432 

U6 

DIP8 

1458 

U7 

DIP8 

LM311 

U8 

DIP8 

LM311 

U9 

DATEL 

TWR-5/1000- 

15/200-D12 

Q1 

TO-92A 

2N2222 
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2. 


Net  List 


NET1 

U4-10 

D902B-A 

NET19 

R411G-2 

R902-1 

R803-2 

R411-1 

NET7  (+15V) 

R803G-2 

U4-11 

U6-8 

NET  12 

U8-3 

R603H-1 

C902-1 

C803-2 

NET2 

C904-1 

C901-1 

C603-2 

R807H-2 

U9-1 

NET20 

U6-3 

R803-1 

D901A-K 

Cl  07-2 

R603-2 

R313-1 

D901B-K 

D107-A 

Ql-3 

C904D-2 

R113C-2 

C808-1 

NET  13  (+5V) 

Ul-7 

NET3 

C708-2 

C907-1 

U6-1 

C306-2 

C514-1 

NET21 

U6-2 

C104-1 

U5-14 

Cl  07-1 

R606-2 

C204-1 

U9-7 

D106-A 

C608-1 

J25-4 

R106-1 

NET4  (-15V) 

C414-1 

R907-2 

Ul-6 

U6-4 

U3-6 

R504-2 

C905-2 

U7-8 

C501-2 

NET22 

C904D-1 

U8-8 

C907G-1 

D107-K 

C305-2 

U2-4 

D106-K 

C804-2 

U4-14 

NET14 

C704-2 

U9-4 

R112-1 

NET23 

Cl  11-2 

Ul-4 

R106-2 

R114-1 

C211-2 

R301-1 

R101-2 

R113A-1 

C408-1 

C807-1 

Ul-1 

R113C-1 

C604-2 

NET8 

R113B-1 

U3-5 

R603-1 

NET  15 

Ul-13 

U7-4 

R603V-2 

R112-2 

U8-4 

R112G-1  . 

NET24 

U2-11 

NET9 

Ul-12 

R114-2 

U4-8 

R603H-2 

Ul-14 

U9-5 

R603V-1 

NET  16 

R311-1 

Ul-11 

R807Q-1 

R704-1 

NET10  (io/10) 

U8-7 

NET25 

R316A-2 

R41 1-2 

R807-2 

R108-1 

R318-2 

R109-1 

R113B-2 

C316-2 

R102A-1 

NET  17 

Ul-8 

U2-13 

R807Q-2 

NET5 

U2-14 

Ql-2 

NET26 

U6-6 

R316B-2 

R108-2 

R607-2 

NET  18 

R 109-2 

R606-1 

NET11 

R807H-1 

Ul-9 

NET6  (Vref) 

C902-2 

U3-14 

U6-7 

C901-2 

R807-1 

NET27 

R607-1 

U9-2 

C807-2 

R110G-1 

R102B-1 

D902A-A 

R110-2 
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Ul-10 

U2-12 

NET28 

NET36 

R102B-2 

C205-1 

R102C-2 

R205-2 

R101-1 

U2-5 

R102A-2 

Ul-2 

NET37 

NET29 

R024J-1 

J25-24 

R113A-2 

U4-7 

NET38 

U4-12 

R203-1 

J25-9 

NET30 

R313-2 

NET39 

U3-13 

R205-1 

R313G-1 

J25-8 

C313-1 

NET31  (vo/50) 

NET40 

R210-1 

R102C-1 

J25-6 

U2-1 

U2-2 

NET41 

NET32  (iL/10) 

R212-1 

J25-7 

R110-1 

U3-12 

NET42 

U2-8 

U3-1 

U2-9 

R301-2 

NET33 

C301-2 

C203-2 

D301-K 

R203-2 

U2-3 

NET43 

NET34 

U3-2 

U3-8 

C210-2 

U3-10 

R210-2 

C302-2 

U2-10 

NET35 

NET44 

U3-4 

C212-2 

J25-13 

R212-2 

NET45 

U3-9 

NET53 

U3-11 

U5-1 

R311-2 

J25-3 

C311-2 

R501-1 

D311-K 

C501-1 

NET46 

NET54 

U3-16 

U5-2 

R316B-1 

J25-2 

R316A-1 

R502-2 

C317-1 

D317-K 

NET55 

U5-3 

NET47 

U5-13 

U3-17 

J25-1 

U7-3 

C317-2 

NET56 

D317-A 

U5-5 

J25-14 

NET48 

R505-2 

U3-18 

R318-1 

NET57 

U5-6 

NET49 

U5-12 

U7-2 

J25-15 

R704-2 

R702-2 

NET58 

C702-2 

U5-9 

U5-10 

NET50 

R509-1 

U7-7 

U5-4 

NET59 

R504-1 

D901A-A 

D902A-K 

NET51 

U8-2 

NET60 

U5-11 

D901B-A 

D902B-K 

NET52  (vin/50) 

U2-6 

NET61 

U2-7 

R907-1 

U4-1 

D907-2 
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PCB  Component  Overlay 


B. 


PHASE-LOCKED  LOOP  PCB  DATA 


This  section  is  divided  into  a  component  parts  listing,  an  EASYTRAX  PCB  netlist,  and  a  PCB 
component  overlay. 


1.  Component  List 


C202 

RAD0.2 

C203 

RAD0.2 

C102 

RAD0.2 

C204 

RAD0.2 

R202 

AXIAL0.4 

C104 

RAD0.2 

R102 

AXIAL0.4 

C108 

RAD0.3 

R107 

AXIAL0.4 

C308 

RAD0.3 

R207 

AXIAL0.4 

C208 

RAD0.3 

C103 

RAD0.2 

JUMPER 

RAD0.2 

2.  Net  List 

NET1 

TSLAVE-2 

C202-1 

R202-2 

NET4 

U2-LM311-2 

R102-1 

TMASTER-2 

NET2 

Cl  02-1 

NET5 

R102-2 

R107-1 

U1-LM31 1-2 

R207-1 

Cl  08-2 

NET3 

C314-2 

R202-1 

C208-2 

R311X 

U1-LM311 

AXIAL0.4 

DIP8 

R311 

U3-4011 

AXIAL0.4 

DIP14 

T-15 

R103 

RAD0.2 

T+15 

AXIAL0.4 

RAD0.2 

TI  SET 

R203 

AXIAL0.4 

RAD0.2 

TSLAVE 

R108 

AXIAL0.4 

RAD0.2 

TMASTER 

C311 

RAD0.2 

RAD0.2 

U2-LM3 1 1 

DIP8 

R208 

AXIAL0.4 

/ 

T+15-2 

NET7 

U2-LM311-8 

R207-2 

U1-LM311-8 

U2-LM31 1-7 

U3-4011-14 

U3-4011-8 

R108-2 

NET8 

R208-1 

C103-1 

U1-LM311-3 

NET6 

R103-2 

R107-2 

R108-1 

U1-LM311-7 

U3-4011-13 

NET9 

C203-1 
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U2-LM311-3 

U1-LM311-4 

NET13 

NET15 

R203-2 

NET11 

R311X-1 

R311-2 

R208-2 

JUMPER-2 

TI_SET-1 

U3-4011-9 

U3-4011-10 

U3-401 1-11 

NET10 

NET14 

C204-2 

NET12 

R311X-2 

Cl  04-2 

JUMPER-1 

R311-1 

T-15-1 

U3-4011-12 

C311-1 

U2-LM311-4 

3.  PCB  Component  Overlay 
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c. 


IGBT  DRIVER  CIRCUIT  PCB  DATA 


This  section  is  divided  into  a  component  parts  listing,  an  EASYTRAX  PCB  netlist,  and  a  PCB 
component  overlay. 


1.  Component  List 


C104,C108 

D123 

2200uF-50V 

1N4148  or  1N914 

RB.3/.6 

DIODEO.4 

C148 

DGATE2,DGATE1 

0.22uF 

1N4744-Zener 

RAD0.3 

DIODEO.4 

C203 

JUMP 

O.luF 

-jumper  with  opto 

RAD0.2 

-input  if  no  opto 

RAD0.1 

D108A,D104A, 

D108B,D104B 

R106 

1N4005 

6.2  -  1/4W 

DIODEO.4 

AXIAL0.4 

R107 

U1 

lk-l/4W 

TLP250  (Toshiba) 

AXIAL0.4 

DIP8 

R201 

U2 

-jumper  with  opto 

HFBR-2521  (HP) 

-1.3k-l/4W  if  no 

HFBR  15/25XX 

opto 

AXIAL0.4 

VR1 

LM78L05ACZ 

R213 

TO-92A 

lk-l/4W 

AXIAL0.4 

2.  Net  List 

NET1 

R201-2 

JUMP-1 

Ul-3 

NET4 

D123-A 

C 148-2 

Cl  04-2 

NET2 

D104A-A 

R213-1 

D104B-A 

U2-3 

Ul-5 

VR1-2 

C203-1 

NET5 

C148-1 

NET3 

C108-1 

R213-2 

D108A-K 

U2-4 

D108B-K 

U2-1 

Ul-8 

VR1-1 

NET9 

Ul-2 

NET6 

R201-1 

D108A-A 

D123-K 

D104A-K 

NET10 

NET7 

R106-1 

D108B-A 

DGATE2-A 

D104B-K 

R107-2 

NET8 

NET11 

Ul-7 

DGATE2-K 

Ul-6 

DGATE1-K 

R106-2 
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PCB  Component  Overlay 


Figure  D-3,  IGBT  Driver  PCB  Component  Overlay 


D.  SENSOR  BOARD  PCB  DATA 

This  section  is  divided  into  a  component  parts  listing,  an  EASYTRAX  PCB 
netlist,  and  a  PCB  component  overlay. 


1.  Component  List 


R401 

J4 

AXIAL0.4 

DIP6 

R201 

J3 

AXIAL0.4 

DIP6 

R301 

U8 

AXIAL0.4 

CL50 

CL50 

R101 

AXIAL0.4 

U7 

C602 

CL50 

CL50 

RAD0.2 

C601 

U6 

CL50 

RAD0.2 

CL50 

C502 

U5 

RAD0.2 

CL50 

CL50 

C501 

RAD0.2 

VR1 

C801 

TO-220 

TI7815 

RAD0.2 

C802 

VR2 

TO-220 

RAD0.2 

MC7915 

C701 

DVR22B 

RAD0.2 

DIODEO.4 

C702 

DVR11B 

RAD0.2 

DIODEO.4 

DVR22A 

R703 

DIODEO.4 

AXIAL0.4 

DVR11A 

R503 

DIODEO.4 

AXIAL0.4 

CVR11 

R803 

RB.3/.6 

AXIAL0.4 

CVR22 

R603 

RB.3/.6 

AXIAL0.4 

U1 

C144 

AD215 

RAD0.2 

AD215 

C344 

U3 

RAD0.2 

AD215 

AD215 

C244 

RAD0.2 

U2 

AD215 

C444 

AD215 

RAD0.2 

U4 

C242 

AD215 

RAD0.2 

AD215 

C442 

J5 

RAD0.2 

DIP6 

C142 

J6 

RAD0.2 

DIP6 

C342 

RAD0.2 
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2. 


Net  List 


NET1 

U6-2 

R401-1 

U5-2 

J4-4 

VR2-3 

Ul-44 

NET2 

U3-44 

R401-2 

U2-44 

U4-1 

U4-44 

C144-1 

NET3 

C344-1 

R201-1 

C244-1 

J4-2 

C444-1 

NET4 

NET10  (+15  V) 

R201-2 

C601-1 

U2-1 

C501-1 

C801-1 

NET5 

C701-1 

R301-1 

U8-1 

U3-1 

U7-1 

U6-1 

NET6 

U5-1 

R301-2 

VR1-3 

J3-4 

Ul-42 

U3-42 

NET7 

U2-42 

R101-1 

U4-42 

J3-1 

C242-1 

C442-1 

NET8 

C142-1 

R101-2 

C342-1 

Ul-1 

NET11 

NET9  (-15  V) 

J4-3 

C602-1 

U2-2 

C502-1 

C802-1 

NET12 

C702-1 

J4-5 

U8-2 

U4-2 

U7-2 

NET13 

J3-2 

NET21 

Ul-2 

DVR22B-K 

DVR11B-A 

NET14 

J3-5 

NET22 

U3-2 

DVR22A-K 

DVR11A-A 

NET  15 

U8-3 

NET23 

J  6-6 

Ul-38 

R803-1 

Jl-3 

NET16 

NET24 

U7-3 

Ul-4 

J5-6 

Ul-3 

R703-1 

NET25 

NET  17 

U3-38 

U6-3 

Jl-4 

J6-1 

R603-1 

NET26 

U3-4 

NET  18 

U3-3 

U5-3 

J5-3 

NET27 

R503-1 

U2-38 

J2-3 

NET  19 

VR1-1 

NET28 

DVR11B-K 

U2-4 

DVR11A-K 

U2-3 

CVR11-1 

NET29 

NET20 

U4-38 

VR2-2 

J2-4 

DYR22B-A 

NET30 

DVR22A-A 

U4-4 

CVR22-2 

U4-3 

147 


PCB  Component  Overlay 


LIST  OF  REFERENCES 


1 .  Blalock,  H.C.,  “Power  Electronic  Converter  Simulation,  Real  Time  Control  and 
Hardware-In-The-Loop  Testing  for  a  Shipboard  Integrated  Power  System,”  Electrical 
Engineer  Thesis,  Naval  Postgraduate  School,  Monterey,  CA,  March,  1995. 

2.  Mohan,  N.,  Undeland,  T.M.,  Robbins,  W.P.,  Power  Electronics,  Converters, 
Applications  and  Design,  John  Wiley  &  Sons,  Inc.,  New  York,  NY,  1989. 

3.  Mayer,  J.  and  Salberta,  F.,  “High-Frequency  Power  Electronic  Converter  for 
Propulsion  Applications,”  Electrical  Engineer  Thesis,  Penn  State  University, 
University  Park,  PA,  January,  1997. 

4.  Arnold,  PC-104K  MPP  Cores,  The  Arnold  Engineering  Company,  Marengo,  IL, 
September,  1995. 

5.  Kassakian,  J.G.,  Schlecht,  M.F.,  and  Verghese,  G.C.,  Principles  of  Power 
Electronics,  Addison-Wesley  Publishing  Company,  New  York,  NY,  1991. 

6.  Ciezki,  J.G.,  “EC2320  Linear  Control  Systems  Class  Notes,”  Naval  Postgraduate 
School,  Monterey,  CA,  October,  1997. 

7.  Hanson,  R.  J.,  “Implementing  Closed-Loop  Control  Algorithms  for  DC-to-DC 
Converters  and  ARCP  Inverters  Using  the  Universal  Controller,”  Electrical  Engineer 
Thesis,  Naval  Postgraduate  School,  Monterey,  CA,  June,  1997. 

8.  Oberley,  M.J.,  “The  Operation  and  Interaction  of  the  Auxiliary  Resonant  Commutated 
Pole  Converter  in  a  Shipboard  DC  Power  Distribution  Network,”  Electrical  Engineer 
Thesis,  Naval  Postgraduate  School,  Monterey,  CA,  December,  1996. 


149 


INITIAL  DISTRIBUTION  LIST 


1 .  Defense  Technical  Information  Center . 2 

8725  John  J.  Kingman  Rd.,  STE  0944 

Fort  Belvoir,  Virginia  22060-6218 

2.  Dudley  Knox  Library . 2 

Naval  Postgraduate  School 

411  Dyer  Rd. 

Monterey,  California  93 943-5 1 0 1 

3.  Chairman,  Code  EC . 1 

Department  of  Electrical  and  Computer  Engineering 

Naval  Postgraduate  School 
Monterey,  California  93943-5121 

4.  Robert  Ashton,  Code  EC/Ah . 3 

Department  of  Electrical  and  Computer  Engineering 

Naval  Postgraduate  School 
Monterey,  California  93943-5121 

5.  John  Ciezki,  Code  EC/Cy . 3 

Department  of  Electrical  and  Computer  Engineering 

Naval  Postgraduate  School 
Monterey,  California  93943-5121 

6.  Roger  Cooley . 2 

Naval  Surface  Warfare  Center 

Carderock  Division/Annapolis  Detachment 
3A  Leggett  Circle 
Annapolis,  MD  21402-5067 

7.  Michael  G.  Badorf. . 2 

c/o  W.P.  Hearn 

R.R.  6  Box  203 
Georgetown,  DE  19947 


